Estimation device, estimation method, and estimation program

ABSTRACT

An estimation device includes at least one processor, in which the processor functions as a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple two-dimensional image acquired by simply imaging the subject. The learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images.

CROSS REFERENCE TO RELATED APPLICATIONS

The present application claims priority under 35 U.S.C. § 119 to Japanese Patent Application No. 2021-039486 filed on Mar. 11, 2021. The above application is hereby expressly incorporated by reference, in its entirety, into the present application.

BACKGROUND Technical Field

The present disclosure relates to an estimation device, an estimation method, and an estimation program.

Related Art

In the related art, energy subtraction processing using two radiation images obtained by irradiating a subject with two types of the radiation having different energy distributions by using an amount of attenuation of transmitted radiation different from each other depending on a substance configuring the subject is known. The energy subtraction processing is a method in which pixels of the two radiation images obtained as described above are associated with each other, and the pixels are multiplied by an appropriate weighting coefficient and then subtracted (subtract) to acquire an image obtained by emphasizing a specific structure. In addition, in addition to a bone part and a soft part, derivation of a composition of a human body, such as fat and muscle, in the soft part is performed by the energy subtraction processing (see JP2018-153605A).

In addition, various methods have been proposed in which the radiation image acquired by imaging the subject is used to derive a radiation image different from the acquired radiation image. For example, U.S. Pat. No. 7,545,965B proposes a method for using a learned model constructed by learning a neural network by using, as teacher data, a radiation image of a subject acquired by simple imaging and a bone part image of the same subject to derive the bone part image from the radiation image of the subject acquired by the simple imaging.

Note that the simple imaging is an imaging method for acquiring one two-dimensional image, which is a transmission image of the subject, by emitting the radiation to the subject once. In the following description, the two-dimensional image acquired by simple imaging will be referred to as a simple two-dimensional image.

However, it is desired to estimate an image in which a specific composition, such as a bone part, is emphasized with higher accuracy.

SUMMARY OF THE INVENTION

The present disclosure has been made in view of the above circumstances, and is to enable estimation of the image in which the specific composition is emphasized with high accuracy.

An aspect of the present disclosure relates to an estimation device comprising at least one processor, in which the processor functions as a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple two-dimensional image acquired by simply imaging the subject, and the learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image for learning may be derived by energy subtraction processing of performing weighting subtraction on the two radiation images.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image may be at least one of a bone part image in which a bone part of the subject is emphasized or a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning may be derived by recognizing the bone part and the soft part of the subject by using at least one radiation image of the two radiation images, deriving attenuation coefficients of the bone part and the soft part by using results of recognition of the bone part and the soft part and the two radiation images, and performing the energy subtraction processing by using the attenuation coefficients.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image may be a bone part image in which a bone part of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning may be derived by deriving new weighting coefficients used for the weighting subtraction based on a pixel value of the bone part included in the bone part image and a pixel value of the soft part included in the soft part image, deriving a new bone part image and a new soft part image by performing the weighting subtraction on the two radiation images by using the new weighting coefficients, and repeating derivation of a further new weighting coefficient based on the new bone part image, derivation of a further new weighting coefficient based on the new soft part image, and derivation of a further new bone part image and a further new soft part image based on the further new weighting coefficients.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image may be a bone part image in which a bone part of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning may be derived by deriving, for each of different energy distributions, a difference between a value of an attenuation coefficient of the soft part×a thickness of the soft part+an attenuation coefficient of the bone part×a thickness of the bone part, and each pixel value of the radiation image while changing, from initial values, the attenuation coefficient of the soft part for each of different energy distributions, the thickness of the soft part, the attenuation coefficient of the bone part for each of different energy distributions, and the thickness of the bone part, deriving the attenuation coefficient of the soft part and the attenuation coefficient of the bone part for each of different energy distributions, at which the difference is minimized or the difference is smaller than a predetermined threshold value, and performing the energy subtraction processing by using a weighting coefficient derived based on the attenuation coefficient of the soft part and the attenuation coefficient of the bone part.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image for learning may be derived by deriving a composition ratio of a plurality of compositions included in a soft part of the subject, deriving, for each pixel of the two radiation images, an attenuation coefficient for each of different energy distributions depending on the composition ratio, and performing the energy subtraction processing by using a weighting coefficient derived based on the derived attenuation coefficient.

In addition, in the estimation device according to the aspect of the present disclosure, the composition ratio may be obtained by deriving a body thickness of the subject as a first body thickness and a second body thickness for each pixel of each of the two radiation images and deriving the composition ratio for each pixel of the radiation image based on the first body thickness and the second body thickness.

In addition, in the estimation device according to the aspect of the present disclosure, the composition ratio may be obtained by deriving the first body thickness and the second body thickness based on an attenuation coefficient of each of the plurality of compositions for each of different energy distributions, deriving the first body thickness and the second body thickness while changing a thickness of the composition and the attenuation coefficient of each composition, and deriving the composition ratio based on the thickness of the composition in which a difference between the first body thickness and the second body thickness is equal to or smaller than a predetermined threshold value.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image for learning may be derived by performing scattered ray removal processing of removing a scattered ray component of the radiation emitted to the subject, which is scattered by the subject, from the two radiation images, and performing the energy subtraction processing on the two radiation images from which the scattered ray component is removed.

In addition, in the estimation device according to the aspect of the present disclosure, the scattered ray removal processing may be performed by acquiring a radiation characteristic of an object interposed between the subject and a radiation detector that detects the radiation image depending on the body thickness distribution, deriving a primary ray distribution and a scattered ray distribution of the radiation included in each of the two radiation images by using the imaging condition, the body thickness distribution, and the radiation characteristic of the object, deriving an error between a sum of the primary ray distribution and the scattered ray distribution of each of the two radiation images and a pixel value at each position of the two radiation images, updating the body thickness distribution such that the error is smaller than a predetermined threshold value, repeating derivation of the radiation characteristic based on the updated body thickness distribution and derivation of the primary ray distribution and the scattered ray distribution included in each of the two radiation images, and subtracting the scattered ray distribution in a case in which the error is smaller than the predetermined threshold value from each of the two radiation images.

In addition, in the estimation device according to the aspect of the present disclosure, the scattered ray removal processing may be performed by deriving a first primary ray distribution and a scattered ray distribution of the radiation transmitted through the subject by using the two radiation images, deriving a second primary ray distribution and a scattered ray distribution of the radiation transmitted through an object interposed between the subject and a radiation detector that detects the radiation image by using the first primary ray distribution, the scattered ray distribution, and a radiation characteristic of the object, and deriving the radiation images after transmission through the subject and the object by using the second primary ray distribution and the scattered ray distribution.

In addition, in the estimation device according to the aspect of the present disclosure, the scattered ray removal processing may be performed by deriving a region detection image by detecting a subject region in which the radiation is transmitted through the subject and reaches a radiation detector and a direct radiation region in which the radiation directly reaches the radiation detector without being transmitted through the subject in the two radiation images, deriving a scattered ray image relating to the scattered ray component based on the region detection image and scattered ray spread information relating to spread of a scattered ray, and subtracting the scattered ray image from the two radiation images.

In addition, in the estimation device according to the aspect of the present disclosure, the emphasis image for learning may be obtained by deriving a processing content of first granulation suppression processing on a first radiation image having higher S/N among the two radiation images, deriving a processing content of second granulation suppression processing on a second radiation image having lower S/N based on the processing content of the first granulation suppression processing, performing granulation suppression processing on the first radiation image based on the processing content of the first granulation suppression processing, performing granulation suppression processing on the second radiation image based on the processing content of the second granulation suppression processing, and deriving the emphasis image for learning by using the two radiation images subjected to the granulation suppression processing.

In addition, in the estimation device according to the aspect of the present disclosure, the processing content of the first granulation suppression processing may be derived based on a physical quantity map of the subject derived based on at least one of the first radiation image or the second radiation image.

Another aspect of the present disclosure relates to an estimation method comprising using a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple radiation image acquired by simply imaging the subject to derive the result of estimation of the at least one emphasis image in which the specific composition of the subject is emphasized from the simple radiation image, in which the learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images.

Note that the estimation method according to the present disclosure may be provided as a program causing a computer to execute.

According to the present disclosure, it is possible to estimate image in which the specific composition is emphasized with high accuracy.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic block diagram showing a configuration of a radiography system to which an estimation device according to a first embodiment of the present disclosure is applied.

FIG. 2 is a diagram showing a schematic configuration of the estimation device according to the first embodiment.

FIG. 3 is a diagram showing a functional configuration of an estimation device according to the first embodiment.

FIG. 4 is a diagram showing a schematic configuration of a neural network used in the present embodiment.

FIG. 5 is a diagram showing teacher data.

FIG. 6 is a diagram showing a schematic configuration of an information derivation device according to the first embodiment.

FIG. 7 is a diagram showing a functional configuration of the information derivation device according to the first embodiment.

FIG. 8 is a diagram showing a bone part image.

FIG. 9 is a diagram showing a soft part image.

FIG. 10 is a diagram for describing learning of the neural network.

FIG. 11 is a conceptual diagram of processing performed by a learned neural network.

FIG. 12 is a diagram showing a display screen of a result of estimation.

FIG. 13 is a flowchart of learning processing performed in the first embodiment.

FIG. 14 is a flowchart showing estimation processing performed in the first embodiment.

FIG. 15 is a diagram showing a functional configuration of an information derivation device according to a second embodiment.

FIG. 16 is a diagram for describing calculation of an index value representing attenuation.

FIG. 17 is a diagram for describing calculation of an index value representing attenuation.

FIG. 18 is a diagram showing a functional configuration of an information derivation device according to a third embodiment.

FIG. 19 is a diagram showing a table that defines a relationship between a body thickness distribution and an initial weighting coefficient.

FIG. 20 is a diagram showing a table that defines a relationship between a pixel value of the bone part and a thickness of the bone part.

FIG. 21 is a diagram showing a table that defines a relationship between a thickness of the soft part, the thickness of the bone part, and a weighting coefficient.

FIG. 22 is a diagram showing a functional configuration of an information derivation device according to a fourth embodiment.

FIG. 23 is a diagram showing a table that defines a relationship between an initial value of the thickness of the soft part and an initial value of an attenuation coefficient of the soft part.

FIG. 24 is a diagram showing a table that defines a relationship between a thickness of the soft part, the thickness of the bone part, and an attenuation coefficient.

FIG. 25 is a diagram showing a functional configuration of an information derivation device according to a fifth embodiment.

FIG. 26 is a diagram for describing a difference between body thicknesses derived from a low-energy image and a high-energy image.

FIG. 27 is a diagram showing a table in which a difference between the body thicknesses derived from two radiation images and a composition ratio of fat are associated with each other.

FIG. 28 is a diagram showing a functional configuration of a scattered ray removal unit in an information derivation device according to a sixth embodiment.

FIG. 29 is a diagram for describing imaging in the sixth embodiment.

FIG. 30 is a diagram for describing the measurement of the scattered ray transmittance depending on the body thickness of the subject.

FIG. 31 is a diagram for describing measurement of a scattered ray transmittance depending on the body thickness of the subject.

FIG. 32 is a table showing a relationship between the body thickness distribution of the subject and the scattered ray transmittance of an object interposed between the subject and a radiation detector.

FIG. 33 is a diagram for describing the measurement of the primary ray transmittance depending on the body thickness of the subject.

FIG. 34 is a diagram for describing measurement of a primary ray transmittance depending on the body thickness of the subject.

FIG. 35 is a table showing a relationship between the body thickness distribution of the subject and the primary ray transmittance of the object interposed between the subject and the radiation detector.

FIG. 36 is a diagram showing a state in which an air layer is interposed between a top plate and a grid.

FIG. 37 is a diagram showing a functional configuration of a scattered ray removal unit in an information derivation device according to a seventh embodiment.

FIG. 38 is a diagram for describing a function of a first derivation unit according to the seventh embodiment.

FIG. 39 is a diagram for describing a method of estimating the body thickness of the subject.

FIG. 40 is a diagram showing a point spread function.

FIG. 41 is a diagram showing a radiation path.

FIG. 42 is a diagram showing a functional configuration of an information derivation device according to an eighth embodiment.

FIG. 43 is a diagram showing a bilateral filter for a first radiation image.

FIG. 44 is a diagram showing a local region of a second radiation image corresponding to a local region of the first radiation image shown in FIG. 43.

FIG. 45 a diagram showing an example of a bilateral filter for a second radiation image.

FIG. 46 is a diagram showing a functional configuration of an information derivation device according to a ninth embodiment.

FIG. 47 is a diagram showing an example of a bilateral filter for a physical quantity map.

FIG. 48 is a diagram showing a functional configuration of a scattered ray removal unit in an information derivation device according to a tenth embodiment.

FIG. 49 is a diagram for describing detection of a subject region and a direct radiation region.

FIG. 50 is a diagram showing a region detection image of a specific line portion.

FIG. 51 is a diagram showing a scattered ray image of the specific line portion.

FIG. 52 is a diagram showing the body thickness distribution of the specific line portion.

FIG. 53 is a diagram showing a functional configuration of an information derivation device according to an eleventh embodiment.

FIG. 54 is a diagram showing an example of energy spectra of the radiation after being transmitted through a muscle tissue and radiation after being transmitted through a fat tissue.

FIG. 55 is a diagram showing still another example of the teacher data.

DETAILED DESCRIPTION

Hereinafter, embodiments of the present disclosure will be described with reference to the drawings. FIG. 1 is a schematic block diagram showing a configuration of a radiography system to which an estimation device according to a first embodiment of the present disclosure is applied. As shown in FIG. 1, the radiography system according to the first embodiment comprises an imaging apparatus 1, an image storage system 9, an estimation device 10 according to the first embodiment, and an information derivation device 50. The imaging apparatus 1, the estimation device 10, and the information derivation device 50 are connected to the image storage system 9 via a network (not shown).

The imaging apparatus 1 is an imaging apparatus that can perform energy subtraction by a so-called one-shot method for converting radiation, such as X-rays, emitted from a radiation source 3 and transmitted through a subject H into energy and irradiating a first radiation detector 5 and a second radiation detector 6 with the converted radiation. At the time of imaging, as shown in FIG. 1, the first radiation detector 5, a radiation energy conversion filter 7 made of a copper plate or the like, and the second radiation detector 6 are disposed in order from a side closest to the radiation source 3, and the radiation source 3 is driven. Note that the first and second radiation detectors 5 and 6 are closely attached to the radiation energy conversion filter 7.

As a result, in the first radiation detector 5, a first radiation image G1 of the subject H by low-energy radiation including so-called soft rays is acquired. In addition, in the second radiation detector 6, a second radiation image G2 of the subject H by high-energy radiation from which the soft rays are removed is acquired. Therefore, the first radiation image G1 and the second radiation image G2 are acquired by imaging the subject H with the radiation having different energy distributions. The first and second radiation images G1 and G2 are input to the estimation device 10. Both the first radiation image G1 and the second radiation image G2 are front images including a periphery of a crotch of the subject H.

The first and second radiation detectors 5 and 6 can perform recording and reading-out of the radiation image repeatedly. A so-called direct-type radiation detector that directly receives emission of the radiation and generates an electric charge may be used, or a so-called indirect-type radiation detector that converts the radiation into visible light and then converts the visible light into an electric charge signal may be used. In addition, as a method for reading out a radiation image signal, it is desirable to use a so-called thin film transistor (TFT) readout method in which the radiation image signal is read out by turning a TFT switch on and off, or a so-called optical readout method in which the radiation image signal is read out by emission of read out light. However, other methods may also be used without being limited to these methods.

In addition, the imaging apparatus 1 can acquire a simple radiation image G0 which is a simple two-dimensional image of the subject H by performing a simple imaging of the subject H by using only the first radiation detector 5. The imaging for acquiring the first and second radiation images G1 and G2 is referred to as energy subtraction imaging in order to distinguish the imaging from simple imaging. In the present embodiment, the first and second radiation images G1 and G2 acquired by the energy subtraction imaging are used as learning data to be described below. In addition, the simple radiation image G0 acquired by the simple imaging is used for deriving the result of estimation of at least one emphasis image in which the specific composition of the subject H is emphasized, as will be described below.

The image storage system 9 is a system that stores image data of the radiation image acquired by the imaging apparatus 1. The image storage system 9 extracts an image corresponding to a request from the estimation device 10 from the stored radiation image and transmits the extracted image to a request source device. Specific examples of the image storage system 9 include picture archiving and communication systems (PACS). Note that in the present embodiment, the image storage system 9 stores a large amount of teacher data for learning the neural network described below.

Then, the estimation device according to the first embodiment will be described. First, a hardware configuration of the estimation device according to the first embodiment will be described with reference to FIG. 2. As shown in FIG. 2, the estimation device 10 is a computer, such as a workstation, a server computer, and a personal computer, and comprises a central processing unit (CPU) 11, a non-volatile storage 13, and a memory 16 as a transitory storage region. In addition, the estimation device 10 comprises a display 14, such as a liquid crystal display, an input device 15, such as a keyboard and a mouse, and a network interface (I/F) 17 connected to a network (not shown). The CPU 11, the storage 13, the display 14, the input device 15, the memory 16, and the network I/F 17 are connected to a bus 18. Note that the CPU 11 is an example of a processor according to the present disclosure.

The storage 13 is realized by a hard disk drive (HDD), a solid state drive (SSD), a flash memory, and the like. The storage 13 as a storage medium stores an estimation program 12A and a learning program 12B installed in the estimation device 10. The CPU 11 reads out the estimation program 12A and the learning program 12B from the storage 13, expands the estimation program 12A and the learning program 12B in the memory 16, and executes the expanded estimation program 12A and the expanded learning program 12B.

Note that the estimation program 12A and the learning program 12B are stored in a storage device of the server computer connected to the network or in a network storage in a state of being accessible from the outside, and are downloaded and installed in the computer that configures the estimation device 10 in response to the request. Alternatively, the estimation program 12A and the learning program 12B are distributed in a state of being recorded on a recording medium, such as a digital versatile disc (DVD) or a compact disc read only memory (CD-ROM), and are installed in the computer that configures the estimation device 10 from the recording medium.

Then, a functional configuration of the estimation device according to the first embodiment will be described. FIG. 3 is a diagram showing the functional configuration of the estimation device according to the first embodiment. As shown in FIG. 3, the estimation device 10 comprises an image acquisition unit 21, an information acquisition unit 22, an estimation unit 23, a learning unit 24, and a display controller 25. Further, the CPU 11 functions as the image acquisition unit 21, the information acquisition unit 22, the estimation unit 23, and the display controller 25 by executing the estimation program 12A. In addition, the CPU 11 functions as the learning unit 24 by executing the learning program 12B.

The image acquisition unit 21 acquires, for example, the first radiation image G1 and the second radiation image G2 which are the front images of the periphery of the crotch of the subject H from the first and second radiation detectors 5 and 6 by causing the imaging apparatus 1 to perform the energy subtraction imaging of the subject H. In a case in which the first radiation image G1 and the second radiation image G2 are acquired, an imaging conditions, such as an imaging dose, a radiation quality, a tube voltage, a source image receptor distance (SID) which is a distance between the radiation source 3 and surfaces of the first and second radiation detectors 5 and 6, a source object distance (SOD) which is a distance between the radiation source 3 and a surface of the subject H., and the presence or absence of a scattered ray removal grid are set.

The SOD and the SID are used to calculate a body thickness distribution as described below. It is preferable that the SOD be acquired by, for example, a time of flight (TOF) camera. It is preferable that the SID be acquired by, for example, a potentiometer, an ultrasound range finder, a laser range finder, or the like.

The imaging conditions need only be set by input from the input device 15 by an operator. The set imaging conditions are stored in the storage 13. The first and second radiation images G1 and G2 acquired by the energy subtraction imaging and the imaging conditions are also transmitted to and stored in the image storage system 9. The first and second radiation images G1 and G2 are used for deriving the teacher data described below.

In addition, the image acquisition unit 21 acquires the simple radiation image G0 which is the front image of the vicinity of the crotch of the subject H by causing the imaging apparatus 1 to perform the simple imaging of the subject H by using only the first radiation detector 5.

Note that in the present embodiment, the first and second radiation images G1 and G2 and the simple radiation image G0 may be acquired by a program separate from the estimation program 12A and stored in the storage 13. In this case, the image acquisition unit 21 acquires the first and second radiation images G1 and G2 stored in the storage 13 by reading out the first and second radiation images G1 and G2 and the simple radiation image G0 from the storage 13 for processing.

The information acquisition unit 22 acquires the teacher data for learning a neural network, which will be described below, from the image storage system 9 via the network I/F 17.

The estimation unit 23 derives the results of estimation of the bone part image in which the bone part is emphasized and the soft part image in which the soft part is emphasized included in the subject H from the simple radiation image G0. Therefore, the estimation unit 23 derives the results of estimation of the bone part image and the soft part image by using a learned neural network 23A that outputs the bone part image and the soft part image in a case in which the simple radiation image G0 is input. Note that in the present embodiment, a target for deriving the results of estimation of the bone part image and the soft part image is an image of the vicinity of the hip joint of the subject H, but the present disclosure is not limited to this.

The learning unit 24 constructs the learned neural network 23A by machine learning the neural network by using the teacher data. Examples of the neural network include a simple perceptron, a multi-layer perceptron, a deep neural network, a convolutional neural network, a deep belief network, a recurrent neural network, and a stochastic neural network. In the present embodiment, the convolutional neural network is used as the neural network.

FIG. 4 is a diagram showing the neural network used in the present embodiment. As shown in FIG. 4, a neural network 30 comprises an input layer 31, an interlayer 32, and an output layer 33. The interlayer 32 comprises, for example, a plurality of convolutional layers 35, a plurality of pooling layers 36, and a fully bonded layer 37. In the neural network 30, the fully bonded layer 37 is present in front of the output layer 33. Further, in the neural network 30, the convolutional layer 35 and the pooling layer 36 are alternately disposed between the input layer 31 and the fully bonded layer 37.

Note that a configuration of the neural network 30 is not limited to the example of FIG. 4. For example, the neural network 30 may comprise one convolutional layer 35 and one pooling layer 36 between the input layer 31 and the fully bonded layer 37.

FIG. 5 is a diagram showing an example of the teacher data used for learning the neural network. As shown in FIG. 5, teacher data 40 consists of learning data 41 and correct answer data 42. In the present embodiment, the data input to the learned neural network 23A in order to obtain the results of estimation of the bone part image and the soft part image is the simple radiation image G0, but the learning data 41 includes two radiation images of the first radiation image G1 and the second radiation image G2 acquired by the energy subtraction imaging.

The correct answer data 42 is the bone part image and the soft part image in the vicinity of the target bone (that is, a femur) of the subject from which the learning data 41 is acquired. The bone part image Gb and the soft part image Gs, which are the correct answer data 42, are derived from the first and second radiation images G1 and G2 by the information derivation device 50. Hereinafter, the information derivation device 50 will be described.

FIG. 6 is a schematic block diagram showing a configuration of the information derivation device according to the first embodiment. As shown in FIG. 6, the information derivation device 50 according to the first embodiment is a computer, such as a workstation, a server computer, and a personal computer, and includes a CPU 51, a non-volatile storage 53, and a memory 56 as a transitory storage region. In addition, the information derivation device 50 includes a display 54, such as a liquid crystal display, an input device 55 including a pointing device, such as a keyboard and a mouse, and a network I/F 57 connected to a network (not shown). The CPU 51, the storage 53, the display 54, the input device 55, the memory 56, and the network I/F 57 are connected to a bus 58.

Like the storage 13, the storage 53 is realized by the HDD, the SSD, the flash memory, and the like. An information derivation program 52 is stored in the storage 53 as the storage medium. The CPU 51 reads out the information derivation program 52 from the storage 53, expands the read out information derivation program 52 in the memory 56, and executes the expanded information derivation program 52.

Then, a functional configuration of the information derivation device according to the first embodiment will be described. FIG. 7 is a diagram showing the functional configuration of the information derivation device according to the first embodiment. As shown in FIG. 7, the information derivation device 50 according to the first embodiment comprises an image acquisition unit 61, a scattered ray removal unit 62, and a subtraction unit 63. Further, the CPU 51 executes the information derivation program 52, so that the CPU 51 functions as the image acquisition unit 61, the scattered ray removal unit 62, and the subtraction unit 63.

The image acquisition unit 61 acquires the first radiation image G1 and the second radiation image G2, which are the learning data 41, stored in the image storage system 9. Note that the image acquisition unit 61 may acquire the first radiation image G1 and the second radiation image G2 by causing the imaging apparatus 1 to image the subject H in the same manner as the image acquisition unit 21 of the estimation device 10.

The image acquisition unit 61 also acquires the imaging conditions in a case in which the first and second radiation images stored in the image storage system 9 are acquired. The imaging conditions include the imaging dose in a case in which the first radiation image G1 and the second radiation image G2 are acquired, the tube voltage, the SID, the SOD, the presence or absence of the scattered ray removal grid, and the like.

Here, each of the first radiation image G1 and the second radiation image G2 includes a scattered ray component based on the radiation scattered in the subject H in addition to a primary ray component of the radiation transmitted through the subject H. Therefore, the scattered ray removal unit 62 removes the scattered ray component from the first radiation image G1 and the second radiation image G2. For example, the scattered ray removal unit 62 may remove the scattered ray component from the first radiation image G1 and the second radiation image G2 by applying a method disclosed in JP2015-043959A. In a case in which a method disclosed in JP2015-043959A or the like is used, the derivation of the body thickness distribution of the subject H and the derivation of the scattered ray component for removing the scattered ray component are performed at the same time.

Hereinafter, the removal of the scattered ray component from the first radiation image G1 will be described, but the removal of the scattered ray component from the second radiation image G2 can also be performed in the same manner. First, the scattered ray removal unit 62 acquires a virtual model of the subject H having an initial body thickness distribution T0(x,y). The virtual model is data virtually representing the subject H of which a body thickness in accordance with the initial body thickness distribution T0(x,y) is associated with a coordinate position of each pixel of the first radiation image G1. Note that the virtual model of the subject H having the initial body thickness distribution T0(x,y) may be stored in the storage 53 of the information derivation device 50 in advance. In addition, the scattered ray removal unit 62 may calculate a body thickness distribution T(x,y) of the subject H based on the SID and the SOD included in the imaging conditions. In this case, the initial body thickness distribution T0(x,y) can be obtained by subtracting the SOD from the SID.

Next, the scattered ray removal unit 62 generates, based on the virtual model, an image obtained by synthesizing an estimated primary ray image in which a primary ray image obtained by imaging the virtual model is estimated and an estimated scattered ray image in which a scattered ray image obtained by imaging the virtual model is estimated as an estimated image in which the first radiation image G1 obtained by imaging the subject H is estimated.

Next, the scattered ray removal unit 62 corrects the initial body thickness distribution T0(x,y) of the virtual model such that a difference between the estimated image and the first radiation image G1 is small. The scattered ray removal unit 62 repeatedly performs the generation of the estimated image and the correction of the body thickness distribution until the difference between the estimated image and the first radiation image G1 satisfies a predetermined termination condition. The scattered ray removal unit 62 derives the body thickness distribution in a case in which the termination condition is satisfied as the body thickness distribution T(x,y) of the subject H. In addition, the scattered ray removal unit 62 removes the scattered ray component included in the first radiation image G1 by subtracting the scattered ray component in a case in which the termination condition is satisfied from the first radiation image G1.

The subtraction unit 63 performs energy subtraction processing to derive a bone part image Gb in which a bone part of the subject H is extracted and a soft part image Gs in which a soft part is extracted from the first and second radiation images G1 and G2. The bone part image Gb and the soft part image Gs derived from the first and second radiation images G1 and G2 by the subtraction unit 63 are examples of an emphasis image for learning in the present disclosure. Note that in the first and second radiation images G1 and G2 in the subsequent processing, the scattered ray component is removed. In a case in which the bone part image Gb is derived, the subtraction unit 63 performs weighting subtraction between the corresponding pixels with respect to the first and second radiation images G1 and G2 as shown in Expression (1) to generate the bone part image Gb in which the bone part of the subject H included in each of the radiation images G1 and G2 is extracted, as shown in FIG. 8. In Expression (1), α is a weighting coefficient. Note that a pixel value of each pixel in a bone region in the bone part image Gb is a bone part pixel value.

Gb(x,y)=α·G2(x,y)−G1(x,y)  (1)

On the other hand, in a case in which the soft part image Gs is derived, the subtraction unit 63 performs weighting subtraction between the corresponding pixels with respect to the first and second radiation images G1 and G2 as shown in Expression (2) to generate the soft part image Gs in which the soft part of the subject H included in each of the radiation images G1 and G2 is extracted, as shown in FIG. 9. In Expression (2), β is a weighting coefficient.

Gs(x,y)=G1(x,y)−β×G2(x,y)  (2)

Note that the soft part image Gs shows a soft region due to a soft tissue of the subject H. In the present embodiment, the “soft tissue” of the subject H refers to a tissue other than a bone tissue, and specifically includes a muscle tissue, a fat tissue, blood, and water.

The bone part image Gb and the soft part image Gs, which are used as the correct answer data 42, are derived at the same time as the time when the learning data 41 is acquired, and are transmitted to the image storage system 9. In the image storage system 9, the learning data 41 and the correct answer data 42 are stored in association with each other as the teacher data 40. Note that in order to improve the robustness of the learning, the teacher data 40 including, as learning data 41, an image obtained by performing at least one of enlargement/reduction, contrast change, movement, in-plane rotation, inversion, or noise addition on the same image may be additionally created and stored.

The description will be returned to the estimation device 10. The learning unit 24 learns the neural network by using a large amount of the teacher data 40. FIG. 10 is a diagram for describing learning of the neural network 30. In a case in which the neural network 30 is learned, the learning unit 24 inputs the learning data 41, that is, the first and second radiation images G1 and G2 to the input layer 31 of the neural network 30. Further, the learning unit 24 outputs the bone part image and the soft part image as output data 47 from the output layer 33 of the neural network 30. Further, the learning unit 24 derives a difference between the output data 47 and the correct answer data 42 as a loss L0. Note that the loss is derived between the bone part image of the output data 47 and the bone part image of the correct answer data 42, and between the soft part image of the output data and the soft part image of the correct answer data 42, respectively, and L0 is used as the reference numeral thereof.

The learning unit 24 learns the neural network 30 based on the loss L0. Specifically, the learning unit 24 adjusts a kernel coefficient in the convolutional layer 35, a weight of the bond between the layers, a weight of the bond in the fully bonded layer 37, and the like (hereinafter referred to as a parameter 48) such that the loss L0 is reduced. For example, an error backpropagation method can be used as a method for adjusting the parameter 48. The learning unit 24 repeats the adjustment of the parameter 48 until the loss L0 is equal to or smaller than a predetermined threshold value. As a result, in a case in which the simple radiation image G0 is input, the parameter 48 is adjusted so as to output the bone part image Gb and the soft part image Gs of the input simple radiation image G0, and the learned neural network 23A is constructed. The constructed learned neural network 23A is stored in the storage 13.

FIG. 11 is a conceptual diagram of processing performed by the learned neural network 23A. As shown in FIG. 11, in a case in which the simple radiation image G0 of a patient is input to the learned neural network 23A constructed as described above, the learned neural network 23A outputs the bone part image Gb and the soft part image Gs for the input simple radiation image G0.

The display controller 25 displays the results of estimation of the bone part image Gb and the soft part image Gs estimated by the estimation unit 23 on the display 14. FIG. 12 is a diagram showing a display screen of the result of estimation. As shown in FIG. 12, the display screen 70 has a first image display region 71 and a second image display region 72. The simple radiation image G0 of the subject H is displayed in the first image display region 71. In addition, in the second image display region 72, the bone part image Gb and the soft part image Gs estimated by the estimation unit 23 are displayed.

Then, the processing performed in the first embodiment will be described. FIG. 13 is a flowchart showing learning processing performed in the first embodiment. First, the information acquisition unit 22 acquires the teacher data 40 from the image storage system 9 (step ST1), and the learning unit 24 inputs the learning data 41 included in the teacher data 40 to the neural network 30 to output the bone part image Gb and the soft part image Gs and learns the neural network 30 by using the loss L0 based on the difference from the correct answer data 42 (step ST2), and the processing returns to step ST1. Further, the learning unit 24 repeats the processing of steps ST1 and ST2 until the loss L0 reaches the predetermined threshold value, and terminates the learning processing. Note that the learning unit 24 may terminate the learning processing by repeating the learning a predetermined number of times. As a result, the learning unit 24 constructs the learned neural network 23A.

Then, estimation processing in the first embodiment will be described. FIG. 14 is a flowchart showing the estimation processing in the first embodiment. Note that the simple radiation image G0 is acquired by the imaging and stored in the storage 13. In a case in which an instruction for starting the processing is input from the input device 15, the image acquisition unit 21 acquires the simple radiation image G0 from the storage 13 (step ST11). Then, the estimation unit 23 derives the results of estimation of the bone part image Gb and the soft part image Gs from the simple radiation image G0 (step ST12). Further, the display controller 25 displays the results of estimation of the bone part image Gb and the soft part image Gs derived by the estimation unit 23 on the display 14 together with the simple radiation image G0 (step ST13), and terminates the processing.

As described above, in the present embodiment, the results of estimation of the bone part image Gb and the soft part image Gs for the simple radiation image G0 are derived by using the learned neural network 23A constructed by performing learning with the first and second radiation images G1 and G2 as teacher data. Here, in the present embodiment, the two radiation images, the first and second radiation images G1 and G2, are used for learning the neural network. Therefore, the learned neural network 23A can derive the results of estimation of the bone part image Gb and the soft part image Gs from the simple radiation image G0 with higher accuracy as compared with a case in which one radiation image and the bone part image Gb and the soft part image Gs are used as the teacher data. Therefore, according to the present embodiment, the results of estimation of the bone part image Gb and the soft part image Gs can be derived with higher accuracy.

In addition, in the first embodiment, in a case in which the bone part image Gb and the soft part image Gs, which are the correct answer data 42, are derived, the bone part image Gb and the soft part image Gs may be derived by recognizing the bone part and the soft part of the subject H by using at least one radiation image of the first radiation image G1 or the second radiation image G2, deriving an attenuation coefficient of the radiation for the bone part and the soft part by using results of recognition of the bone part and the soft part and the first and second radiation images G1 and G2, performing the energy subtraction processing by using the derived attenuation coefficient. Hereinafter, this case will be described as a second embodiment. Note that the energy subtraction processing in the second embodiment is disclosed in, for example, WO2020/175319A.

FIG. 15 is a diagram showing a functional configuration of an information derivation device according to the second embodiment. Note that in FIG. 15, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. As shown in FIG. 15, an information derivation device 50A according to the second embodiment further comprises a structure recognition unit 65 and a weighting coefficient derivation unit 66 with respect to the information derivation device 50 according to the first embodiment.

The structure recognition unit 65 recognizes a structure included in the subject H by using at least one of the first radiation image G1 or second radiation image G2. In the second embodiment, the structure to be recognized is the bone part and the soft part. In the second embodiment, the structure recognition unit 65 uses both the first and second radiation images G1 and G2 acquired by the image acquisition unit 61 for recognition processing, but any one of the first radiation image G1 or the second radiation image G2 may be used for the recognition processing. Note that the structure recognition unit 65 may recognize the structure by using the first and second radiation images G1 and G2 after scattered ray removal processing, or may recognize the structure by using the first and second radiation images G1 and G2 before the scattered ray removal processing.

The structure recognition unit 65 recognizes a position, a size, and/or a shape of the structure included in the subject H appearing in the first and second radiation images G1 and G2. That is, the recognition processing performed by the structure recognition unit 65 is processing of specifying the position, the size, and/or the shape of the structure having a boundary with another tissue or the like in the subject H appearing in the radiation image. In the second embodiment, the bone part and the soft part of the subject H are recognized as the structure.

The weighting coefficient derivation unit 66 derives an index value representing attenuation of the radiation as weighting coefficients α and β to be used for the subtraction processing by using the result of recognition of the structure recognition unit 65 and the first and second radiation images G1 and G2, for the structure recognized by the structure recognition unit 65. The attenuation coefficient is a so-called linear attenuation coefficient, and represents a degree (ratio) of the radiation attenuation due to absorption or scattering. The attenuation coefficient differs depending on a specific composition (density or the like) and the thickness (mass) of the structure through which radiation is transmitted.

The weighting coefficient derivation unit 66 calculates the index value representing the attenuation by using a ratio or a difference of the pixel values in the corresponding pixels of the first and second radiation images G1 and G2. FIGS. 16 and 17 are diagrams for describing the calculation of the index value representing the attenuation. Note that in FIG. 16 and FIG. 17, the calculation of the index value representing the attenuation for the three structures will be described. As shown in FIG. 16, it is assumed that three types of structures having the compositions “Ca”, “Cb”, and “Cc” are present in the first radiation image G1 and these pixel values in the first radiation image G1 are all “V1”. On the other hand, in the second radiation image G2, the pixel values of the corresponding pixels are “Va”, “Vb”, and “Vc”, respectively. A degree of decrease in the pixel value corresponds to a degree of the attenuation of the radiation by each structure (each composition). Therefore, as shown in FIG. 17, the weighting coefficient derivation unit 66 can calculate an index value μa representing the attenuation of the structure of the composition “Ca”, an index value μb representing the attenuation of the structure of the composition “Cb”, and an index value μc representing the attenuation of the structure of the composition “Cc” by using the ratio or the difference of the pixel value of the first radiation image G1 and the pixel value of the corresponding second radiation image G2. In the second embodiment, the weighting coefficient derivation unit 66 calculates the index value representing the attenuation of the bone part and the index value representing the attenuation of the soft part by using the bone part and the soft part as the composition.

Note that in a case in which the ratio or the difference of the pixel value of the first radiation image G1 and the pixel value of the corresponding pixel of the second radiation image G2 is known, the index value μ representing the attenuation can be calculated. Therefore, for the sake of description, the pixel values of the compositions “Ca”, “Cb”, and “Cc” are set to be “V1” common in the first radiation image G1, but it is not necessary to set the pixel values of the compositions “Ca”, “Cb”, and “Cc” to be common in the first radiation image G1.

In the present embodiment, the bone part image Gb and the soft part image Gs are derived by the subtraction unit 63. Therefore, the weighting coefficient derivation unit 66 derives a ratio Gs1(x,y)/Gs2(x,y) of a pixel value Gs1(x,y) of the first radiation image G1, which is the low-energy image, and a pixel value Gs2(x,y) of the second radiation image G2, which is the high-energy image, as the index value representing the attenuation, that is, the weighting coefficient α in Expression (1) for the soft region in the first and second radiation images G1 and G2. Note that the ratio Gs1(x,y)/Gs2(x,y) represents a ratio μls/μhs of an attenuation coefficient μls for the low-energy radiation to an attenuation coefficient μhs for the high-energy radiation in the soft part.

In addition, the weighting coefficient derivation unit 66 derives a ratio Gb1(x,y)/Gb2(x,y) of a pixel value Gb1(x,y) of the first radiation image G1, which is the low-energy image, and a pixel value Gb2(x,y) of the second radiation image G2, which is the high-energy image, as the index value representing the attenuation coefficient, that is, the weighting coefficient β in Expression (2) for the bone region in the first and second radiation images G1 and G2. Note that the ratio Gb1 (x,y)/Gb2 (x,y) represents the ratio μlb/μhb of the attenuation coefficient μlb for the low-energy radiation to the attenuation coefficient μhb for the high-energy radiation of the bone part.

In the second embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by performing the energy subtraction processing by Expression (1) and Expression (2) by using the weighting coefficients α and β derived by the weighting coefficient derivation unit 66. The weighting coefficients α and β in Expression (1) and Expression (2) is derived from the attenuation coefficient of the bone part and the attenuation coefficient of the soft part derived by the weighting coefficient derivation unit 66.

Then, a third embodiment of the present disclosure will be described. FIG. 18 is a diagram showing a functional configuration of an information derivation device according to the third embodiment. Note that in FIG. 18, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. In the third embodiment, in a case in which the bone part image Gb and the soft part image Gs, which are the correct answer data 42, are derived, the bone part image Gb and the soft part image Gs are derived by deriving new weighting coefficients used for the weighting subtraction based on the pixel value of the bone part included in the bone part image Gb and the pixel value of the soft part included in the soft part image Gs, deriving a new bone part image and a new soft part image by performing the weighting subtraction on the first and second radiation images G1 and G2 by using the new weighting coefficients, and repeating the derivation of a further new weighting coefficient based on the new bone part image and the new soft part image and the derivation of a further new bone part image and a further new soft part image based on the further new weighting coefficients.

Therefore, as shown in FIG. 18, an information derivation device 50B according to the third embodiment further comprises an initial weighting coefficient setting unit 67 and a weighting coefficient derivation unit 68 with respect to the information derivation device 50 according to the first embodiment.

The initial weighting coefficient setting unit 67 sets the initial weighting coefficient, which is the initial value of the weighting coefficient in a case in which the subtraction unit 63 performs the subtraction processing, based on the body thickness distribution in a case in which the scattered ray removal unit 62 satisfies the termination condition. Here, in the third embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by Expression (1) and Expression (2) by using the initial weighting coefficient set by the initial weighting coefficient setting unit 67 and the weighting coefficients α and β derived by the weighting coefficient derivation unit 68.

The initial weighting coefficient setting unit 67 sets initial weighting coefficients α0 and β0, which are the initial value of the weighting coefficient, based on the body thickness distribution in a case in which the scattered ray removal unit 62 satisfies the termination condition. In the third embodiment, as shown in FIG. 19, a table LUT2 that defines a relationship between the body thickness distribution and the initial weighting coefficients α0 and β0 is stored in the storage 53. The initial weighting coefficient setting unit 67 sets the initial weighting coefficients α0 and β0 based on the body thickness distribution with reference to the table LUT2.

Here, since a degree of the beam hardening described above depends on a thickness ts of the soft part and a thickness tb of the bone part in the subject H, an attenuation coefficient μs of the soft part and an attenuation coefficient μb of the bone part can be defined as μs(ts, tb) and μb(ts, tb) as functions of ts and tb.

In the energy subtraction processing, since there are the images having two different energy distributions, the attenuation coefficient of the soft part of the low-energy image (first radiation image G1 in the present embodiment) can be represented by μls(ts, tb), and the attenuation coefficient of the bone part thereof can be represented by μlb(ts, tb). In addition, the attenuation coefficient of the soft part of the high-energy image (second radiation image G2 in the present embodiment) can be represented by μhs(ts, tb), and the attenuation coefficient of the bone part thereof can be represented by μhb(ts, tb).

In order to derive the bone part image Gb, it is necessary to eliminate the contrast of the soft part included in the radiation image. Therefore, the weighting coefficient α can be obtained by α=μls(ts, tb)/μhs(ts, tb) by using the ratio of the attenuation coefficient of the soft part. In addition, in order to derive the soft part image Gs, it is necessary to eliminate the contrast of the bone part included in the radiation image. Therefore, the weighting coefficient β can be obtained by β=μlb(ts, tb)/μhb(ts, tb) using the ratio of the attenuation coefficient of the bone part. That is, the weighting coefficients α and β can be represented as a function of the thickness ts of the soft part and the thickness tb of the bone part.

In the third embodiment, the subtraction unit 63 first performs the subtraction processing of performing the weighting subtraction on the first and second radiation images G1 and G2 between the corresponding pixels by using the initial weighting coefficients α0 and β0 set by the initial weighting coefficient setting unit 67. Thereafter, as will be described below, the subtraction processing is performed by using weighting coefficients αnew and βnew derived by the weighting coefficient derivation unit 68.

The weighting coefficient derivation unit 68 derives the new weighting coefficients αnew and βnew based on the pixel value Gb(x,y) of the bone part included in the bone part image Gb. Here, the pixel value Gb(x,y) of the bone part corresponds to the thickness of the bone part of the subject H. Therefore, in the third embodiment, the radiation image of a standard object is acquired as a standard radiation image by imaging the standard object simulating the bone part having various thicknesses in advance. Further, in advance, a table that defines the relationship between the pixel value of the bone part and the thickness thereof is derived by using a relationship between the pixel value of the region of the standard object and the thickness of the standard object in the standard radiation image and stored in the storage 53. FIG. 20 is a diagram showing the table that defines the relationship between the pixel value of the bone part and the thickness of the bone part. In a table LUT3 shown in FIG. 20, the bone part is thicker as the pixel value Gb(x,y) of the bone part is lower (that is, the brightness is higher).

The weighting coefficient derivation unit 68 derives the thickness tb of the bone part in each pixel of the bone part image Gb from each pixel value Gb(x,y) of the bone part image Gb with reference to the table LUT3. Note that the region in the bone part image Gb in which the bone part is not present consists only of the soft part, the thickness tb of the bone part is 0. On the other hand, in the bone part image Gb, in the pixel in which the thickness tb of the bone part is not 0, the weighting coefficient derivation unit 68 derives the thickness ts of the soft part by subtracting the thickness tb of the bone part from the body thickness distribution in a case in which the scattered ray removal unit 62 satisfies the termination condition.

As described above, the weighting coefficients α and β can be represented as the function of the thickness ts of the soft part and the thickness tb of the bone part. In the third embodiment, the storage 53 stores a table that defines a relationship between the thickness ts of the soft part, the thickness tb of the bone part, and the weighting coefficients α and β. FIG. 21 is a diagram showing a table that defines the relationship between the thickness ts of the soft part, the thickness tb of the bone part, and the weighting coefficients α and β. As shown in FIG. 21, a table LUT4 three-dimensionally represents the relationship between the thickness ts of the soft part, the thickness tb of the bone part, and the weighting coefficient α (or β). Here, in the table LUT4, the weighting coefficient α (or β) is smaller as the thickness ts of the soft part and the thickness tb of the bone part are larger.

Note that in the third embodiment, a plurality of the tables LUT4 are prepared and stored in the storage of the information derivation device 50B in accordance with the energy distribution of the radiation used at the time of imaging. The weighting coefficient derivation unit 68 acquires information on the energy distribution of the radiation used at the time of imaging based on the imaging conditions, reads out the table LUT4 corresponding to the acquired information of the energy distribution from the storage, and uses the read out table LUT4 for deriving the weighting coefficient. Further, the weighting coefficient derivation unit 68 derives the new weighting coefficients αnew and βnew with reference to the table LUT4 based on the derived thickness tb of the bone part and the thickness ts of the soft part.

The subtraction unit 63 derives a new bone part image Gbnew and a new soft part image Gsnew by Expression (1) and Expression (2) by using the new weighting coefficients αnew and βnew derived by the weighting coefficient derivation unit 68.

Note that the new bone part image Gbnew and the new soft part image Gsnew may be used as the final bone part image Gb and soft part image Gs, but in the third embodiment, the derivation of the weighting coefficients α and β and the subtraction processing are repeatedly performed.

That is, the weighting coefficient derivation unit 68 derives a new thickness tbnew of the bone part with reference to the table LUT3 based on the pixel value of the bone part in the new bone part image Gbnew. Further, the weighting coefficient derivation unit 68 derives a difference Δtb between a new thickness tbnew of the bone part and the thickness tb of the bone part obtained in the previous processing, and determines whether or not the difference Δtb is smaller than a predetermined threshold value. In a case in which the difference Δtb is equal to or larger than the threshold value, the weighting coefficient derivation unit 68 derives a new thickness tsnew of the soft part from the new thickness tbnew of the bone part, and derives further new weighting coefficients αnew and βnew with reference to the table LUT4 based on the new thickness tbnew of the bone part and the new thickness tsnew of the soft part.

The subtraction unit 63 performs the subtraction processing by using the further new weighting coefficients αnew and βnew, and derives a further new bone part image Gbnew and a further new soft part image Gsnew.

Then, the weighting coefficient derivation unit 68 derives a further new thickness tbnew of the bone part based on the further new bone part image Gbnew, and derives a difference Δtb between the further new thickness tbnew of the bone part and the thickness tb of the bone part obtained in the previous processing.

In the third embodiment, the subtraction unit 63 and the weighting coefficient derivation unit 68 repeat the subtraction processing and the derivation of the weighting coefficients αnew and βnew until the difference Δtb derived by the weighting coefficient derivation unit 68 is smaller than the predetermined threshold value.

Further, in the third embodiment, the bone part image Gb and the soft part image Gs in a case in which the difference Δtb is smaller than the threshold value are used as the correct answer data.

Then, a fourth embodiment of the present disclosure will be described. FIG. 22 is a diagram showing a functional configuration of an information derivation device according to the fourth embodiment. Note that in FIG. 22, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. In the fourth embodiment of the present disclosure, in a case in which the bone part image Gb and the soft part image Gs, which are the correct answer data 42, are derived, the bone part image Gb and the soft part image Gs are derived by deriving, for each of different energy distributions, the difference between a value of the attenuation coefficient of the soft part×the thickness of the soft part+the attenuation coefficient of the bone part×the thickness of the bone part, and each pixel value of the radiation image while changing, from initial values, the attenuation coefficient of the soft part for each of different energy distributions, the thickness of the soft part, the attenuation coefficient of the bone part for each of different energy distributions, and the thickness of the bone part, deriving the attenuation coefficient of the soft part and the attenuation coefficient of the bone part for each of different energy distributions, at which the difference is minimized or the difference is smaller than the predetermined threshold value, and performing the energy subtraction processing by using the weighting coefficient derived based on the attenuation coefficient of the soft part and the attenuation coefficient of the bone part.

Therefore, as shown in FIG. 22, an information derivation device 50C according to the fourth embodiment further comprises an initial value derivation unit 81, an attenuation coefficient derivation unit 82, and a weighting coefficient derivation unit 83 with respect to the information derivation device 50 according to the first embodiment.

The initial value derivation unit 81 derives the initial values of the attenuation coefficient of the soft part for each of different energy distributions, the thickness of the soft part, the attenuation coefficient of the bone part for each of different energy distributions, and the thickness of the bone part for deriving the weighting coefficient in a case of performing the energy subtraction processing. Specifically, the initial values μls0, μhs0, ts0, μlb0, μhb0, and tb0 of the attenuation coefficient μls of the soft part for the low-energy radiation, the attenuation coefficient μhs of the soft part for the high-energy radiation, the thickness ts of the soft part, the attenuation coefficient μlb of the bone part for the low-energy radiation, the attenuation coefficient μhb of the bone part for the high-energy radiation, and the thickness tb of the bone part are derived.

In the fourth embodiment, as shown in Expression (1) and Expression (2), the subtraction unit 63 derives the bone part image Gb obtained by extracting the bone part in the subject H and the soft part image Gs obtained by extracting the soft part by performing the subtraction processing of performing the weighting subtraction on the first and second radiation images G1 and G2 between the corresponding pixels by using the weighting coefficient derived by the weighting coefficient derivation unit 83 as described below.

Here, as described above, the attenuation coefficient μs of the soft part and the attenuation coefficient μb of the bone part can be defined as μs(ts, tb) and μb(ts, tb) as functions of ts and tb. Therefore, in the fourth embodiment as well as in the third embodiment, the attenuation coefficient of the soft part of the low-energy image (first radiation image G1 in the present embodiment) can be represented by μls(ts, tb), and the attenuation coefficient of the bone part can be represented by μlb(ts, tb). In addition, the attenuation coefficient of the soft part of the high-energy image (second radiation image G2 in the present embodiment) can be represented by μhs(ts, tb), and the attenuation coefficient of the bone part thereof can be represented by μhb(ts, tb).

In order to derive the bone part image Gb, it is necessary to eliminate the contrast of the soft part included in the radiation image. Therefore, the weighting coefficient α can be obtained by α=μls(ts, tb)/μhs(ts, tb) by using the ratio of the attenuation coefficient of the soft part. In addition, in order to derive the soft part image Gs, it is necessary to eliminate the contrast of the bone part included in the radiation image. Therefore, the weighting coefficient β can be obtained by β=μlb(ts, tb)/μhb(ts, tb) using the ratio of the attenuation coefficient of the bone part. Note that in the following description, the attenuation coefficients μls(ts, tb), μhs(ts, tb), μlb(ts, tb), and μhb(ts, tb) will be simply represented by the attenuation coefficients μls, μhs, μlb, and μhb by omitting (ts, tb).

The initial value derivation unit 81 uses the body thickness distribution in a case in which the scattered ray removal unit 62 satisfies the termination condition as the initial value ts0 of the thickness ts of the soft part. Note that in the fourth embodiment, the body thickness distribution used in a case of performing the scattered ray removal processing is the body thickness distribution assuming that the subject H consists of only the soft part. Therefore, the initial value tb0 of the thickness tb of the bone part is 0. In addition, as the initial values μls0, μhs0, μlb0, and μhb0 of the attenuation coefficients, values corresponding to the initial values ts0 and tb0 of the thicknesses ts and tb of the soft part and bone part are derived. In the fourth embodiment, since the initial value tb0 of the thickness tb of the bone part is 0, the attenuation coefficients μlb0 and μhb0 of the bone part are 0. As the attenuation coefficients μls0 and μhs0 of the soft part, a value corresponding to the initial value ts0 of the thickness ts of the soft part is derived. Therefore, in the fourth embodiment, a table that defines a relationship between the initial value ts0 of the thickness ts of the soft part and the initial values μls0 and μhs0 of the attenuation coefficients of the soft part is stored in the storage 53.

FIG. 23 is a diagram showing a table that defines the relationship between the initial value ts0 of the thickness ts of the soft part and the initial values μls0 and μhs0 of the attenuation coefficients of the soft part. The initial value derivation unit 81 derives the initial values μls0 and μhs0 of the attenuation coefficient of the soft part depending on the initial value ts0 of the thickness ts of the soft part with reference to a table LUT5 stored in the storage 53.

The attenuation coefficient derivation unit 82 derives the attenuation coefficients μls and μhs of the soft part and the attenuation coefficients μlb and μhb of the bone part for each of different energy distributions. Here, for the energy subtraction processing, the low-energy image and the high-energy image are acquired by imaging the subject H with the radiation having different energy distributions. In the present embodiment, the first radiation image G1 is the low-energy image and the second radiation image G2 is the high-energy image. The pixel value G1(x,y) of each pixel of the first radiation image G1, which is the low-energy image, and the pixel value G2(x,y) of each pixel of the second radiation image G2, which is the high-energy image, are represented by Expression (3) and Expression (4) by using the thickness ts(x,y) of the soft part, the thickness tb(x,y) of the bone part, and the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) at the corresponding pixel position. Note that in Expression (3) and Expression (4), the description of (x,y) is omitted.

G1=μls×ts+μlb×tb  (3)

G2=μhs×ts+μhb×tb  (4)

In order to derive the weighting coefficients α and β for performing the energy subtraction processing, it is necessary to derive the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y). The attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) are represented as functions of the thickness ts of the soft part and the thickness tb of the bone part as described above. Therefore, in order to derive the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y), it is necessary to derive the thickness ts of the soft part and the thickness tb of the bone part. ts and tb in Expression (3) and Expression (4) are solved by Expression (5) and Expression (6).

ts={μhb×G1−μlb×G2}/{μls×μhb−μlb×μhs}  (5)

tb={μls×G2−μhs×G1}/{μls×μhb−μlb×μhs}  (6)

Here, the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) on the right side of Expression (5) and Expression (6) are represented as the functions of the thickness ts of the soft part and the thickness tb of the bone part, and thus Expression (5) and Expression (6) cannot be solved algebraically.

Therefore, in the fourth embodiment, error functions EL and EH shown in Expression (7) and Expression (8) are set. The error functions EL and EH correspond to the difference between the value of the attenuation coefficient of the soft part×the thickness of the soft part+the attenuation coefficient of the bone part×the thickness of the bone part and each pixel value of the radiation image for each of different energy distributions. Further, in order to minimize the error functions EL and EH at the same time, in the fourth embodiment, an error function E0 shown in Expression (9) is set. Further, while changing the thickness ts of the soft part, the thickness tb of the bone part, and the attenuation coefficients μls, μhs, μlb, and μhb from the initial values, the error function E0 is minimized, or a combination of the thickness ts of the soft part and the thickness tb of the bone part at which the error function E0 is smaller than a predetermined threshold value is derived. In this case, it is preferable to derive the thickness ts of the soft part and the thickness tb of the bone part by using an optimization algorithm, such as the steepest descent method and the conjugate gradient method. As the initial values of the thickness ts of the soft part, the thickness tb of the bone part, and the attenuation coefficients μls, μhs, μlb, and μhb used in this case, ts0, tb0, μls0, μhs0, μlb0, and μhb0 derived by the initial value derivation unit 81 are used.

EL=G1−{μls×ts+μlb×tb}  (7)

EH=G2−{μhs×ts+μhb×tb}  (8)

E0=EL ² +Eh ²  (9)

Note that the attenuation coefficient used in the process for deriving the thickness ts of the soft part and the thickness tb of the bone part is derived with reference to a table that defines a relationship between the thickness ts of the soft part, the thickness tb of the bone part, which are predetermined, and the attenuation coefficient. Such a table is stored in the storage 53 of the information derivation device.

FIG. 24 is a diagram showing the table that defines the relationship between the thickness ts of the soft part, the thickness tb of the bone part, and the attenuation coefficient. As shown in FIG. 24, the table LUT6 three-dimensionally represents the relationship between the thickness ts of the soft part and the thickness tb of the bone part and the attenuation coefficient μ. Note that although only one LUT6 is shown in FIG. 24, the table is prepared for each of the attenuation coefficients μls, μhs, μlb, and μhb and stored in the storage. Here, in the table LUT6, the attenuation coefficient μ is smaller as the thickness ts of the soft part and the thickness tb of the bone part are larger.

In a case in which the thickness ts of the soft part and the thickness tb of the bone part in a case in which the error function E0 is minimized or the error function E0 is smaller than the predetermined threshold value are derived, the attenuation coefficient derivation unit 82 derives the attenuation coefficients μls, μhs, μlb, and μhb with reference to the table LUT6.

The weighting coefficient derivation unit 83 derives the weighting coefficients α and β used by the subtraction unit 63 in a case of performing the subtraction processing. That is, the weighting coefficient derivation unit 83 derives the weighting coefficients α and β by performing a calculation of α=μls/μhs and β=μlb/μhb by using the attenuation coefficients μls, μhs, μlb, and μhb derived by the attenuation coefficient derivation unit 82.

In the fourth embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by Expression (1) and Expression (2) by using the weighting coefficients α and β derived by the weighting coefficient derivation unit 83.

Then, a fifth embodiment of the present disclosure will be described. FIG. 25 is a diagram showing a functional configuration of an information derivation device according to the fifth embodiment. Note that in FIG. 25, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. In the fifth embodiment of the present disclosure, in a case in which the bone part image Gb and the soft part image Gs, which are the correct answer data 42, are derived, the bone part image Gb and the soft part image Gs are derived by deriving a composition ratio of a plurality of compositions included in the soft part of the subject H, deriving, for each pixel of the first and second radiation images G1 and G2, the attenuation coefficient of the soft part for each of different energy distributions depending on the composition ratio, and performing the subtraction processing by using the weighting coefficient derived based on the derived attenuation coefficient of the soft part and the predetermined attenuation coefficient of the bone part.

Therefore, as shown in FIG. 25, an information derivation device 50D according to the fifth embodiment further comprises a composition ratio derivation unit 84 and an attenuation coefficient setting unit 85 with respect to the information derivation device 50 according to the first embodiment.

The composition ratio derivation unit 84 acquires the composition ratio of the subject H. In the fifth embodiment, the composition ratio derivation unit 84 acquires the composition ratio by deriving the composition ratio of the subject H based on the first and second radiation images G1 and G2. Note that in the fifth embodiment, the composition ratio of the fat is derived as the composition ratio. Therefore, although the subject H includes the bone part, for the sake of description, the description will be made on assumption that the first and second radiation images G1 and G2 do not include the bone part and include only the soft part.

In order to derive the composition ratio, the composition ratio derivation unit 84 first derives the body thicknesses of the subject H for each pixel of the first and second radiation images G1 and G2 as the first body thickness and the second body thickness. Specifically, the composition ratio derivation unit 84 derives a first body thickness t1 of the subject H by assuming that a brightness distribution of the first radiation image G1 coincides with the body thickness distribution of the subject H, and converting the pixel value of the first radiation image G1 into the thickness by using the attenuation coefficient of the radiation in the muscle of the subject H. In addition, the composition ratio derivation unit 84 derives a second body thickness t2 of the subject H by assuming that a brightness distribution of the second radiation image G2 coincides with the body thickness distribution of the subject H, and converting the pixel value of the second radiation image G2 into the thickness by using the attenuation coefficient in the muscle of the subject H.

Here, since the degree of the beam hardening described above depends on a thickness tf of the fat and a thickness tm of the muscle in the subject H, an attenuation coefficient μf of the fat and an attenuation coefficient μm of the muscle are can be defined as μf(tf, tm) and μm(tf, tm) as non-linear functions of the thickness tf of the fat and the thickness tm of the muscle.

In the fifth embodiment, the attenuation coefficient of the fat of the first radiation image G1, which is the low-energy image, can be represented as μlf(tf, tm), and the attenuation coefficient of the muscle can be represented as μlm(tf, tm). In addition, the attenuation coefficient of the fat of the second radiation image G2, which is the high-energy image, can be represented as μhf(tf, tm), and the attenuation coefficient of the muscle can be represented as μhm(tf, tm).

In addition, the pixel value G1(x,y) of each pixel in the soft region of the first radiation image G1, which is the low-energy image, and the pixel value G2(x,y) of each pixel in the soft region of the second radiation image G2, which is the high-energy image, are represented by Expression (10) and Expression (11), respectively, by using the thickness tf(x,y) of the fat, the thickness tm(x,y) of the muscle, and the attenuation coefficients μlf(x,y), μhf(x,y), μlm(x,y), and μhm(x,y) at the corresponding pixel position. Note that in Expression (10) and Expression (11), the description of (x,y) is omitted.

G1=μlf×tf+μlm×tm  (10)

G2=μhf×tf+μhm×tm  (11)

As described above, in the fifth embodiment, in a case of deriving the first body thickness t1 and the second body thickness t2, the pixel values of the first radiation image G1 and the second radiation image G2 are converted into the thicknesses by using the attenuation coefficient of the muscle in the subject H. Therefore, in the fifth embodiment, the composition ratio derivation unit 84 derives the first body thickness t1 and the second body thickness t2 by Expression (12) and Expression (13). Note that in Expression (12) and Expression (13), the description of (x,y) is omitted.

t1=G1/μlm  (12)

t2=G2/μhm  (13)

In a case in which the subject H includes only the muscle at the pixel positions from which the first and second body thicknesses t1 and t2 are derived, the first body thickness t1 and the second body thickness t2 coincide with each other. However, in the actual subject H, both the muscle and the fat are included at the same pixel positions of the first and second radiation images G1 and G2. Therefore, the first and second body thicknesses t1 and t2 derived by Expression (12) and Expression (13) do not coincide with the actual body thickness of the subject H. In addition, in the first body thickness t1 derived from the first radiation image G1, which is the low-energy image, and the second body thickness t2 derived from the second radiation image G2, which is the high-energy image, the first body thickness t1 has a larger value than the second body thickness t2.

For example, as shown in FIG. 26, it is assumed that the actual body thickness is 100 mm and the thicknesses of the fat and the muscle are 30 mm and 70 mm, respectively. In this case, the first body thickness t1 derived from the first radiation image G1 acquired by the low-energy radiation is derived to be, for example, 80 mm, and the second body thickness t2 derived by the second radiation image G2 acquired by the high-energy radiation is derived to be, for example, 70 mm. In addition, a difference between the first body thickness t1 and the second body thickness t2 is larger as the composition ratio of the fat is larger.

Here, the difference between the first body thickness t1 and the second body thickness t2 is changed depending on the composition ratio of the fat and the muscle in the subject H. Therefore, in the fifth embodiment, a subject model in which the composition ratio of the fat is variously changed is imaged with the radiation having different energy distributions, the body thicknesses are derived from the two radiation images acquired in this way, respectively, and a table in which a difference between the body thicknesses derived from the two radiation images and the composition ratio of the fat are associated with each other is created in advance and stored in the storage 53.

FIG. 27 is a diagram showing the table in which the difference between the body thicknesses derived from the two radiation images and the composition ratio of the fat are associated with each other. As shown in FIG. 27, in a table LUT7, a horizontal axis is the difference between the body thicknesses derived from the two radiation images, and a vertical axis is the composition ratio of the fat. As shown in FIG. 27, the composition ratio of the fat is larger as the larger the difference between the body thicknesses derived from the two radiation images is larger. Note that the table in which the difference between the body thicknesses derived from the two radiation images and the composition ratio of the fat are associated with each other is prepared for each radiation energy distribution used at the time of imaging and stored in the storage 53.

The composition ratio derivation unit 84 derives the difference between the first body thickness t1 and the second body thickness t2, which are derived, and derives a composition ratio R(x,y) of the fat with reference to the LUT7 stored in the storage of the information derivation device 50D. Note that the composition ratio of the muscle can be derived by subtracting the derived composition ratio R(x,y) of the fat from 100%.

The attenuation coefficient setting unit 85 sets the attenuation coefficient of the radiation used in a case of acquiring the first and second radiation images G1 and G2 for each pixel of the first and second radiation images G1 and G2 depending on the composition ratio R(x,y) of the fat. Specifically, the attenuation coefficient of the soft part of the subject H is set. Here, in the fifth embodiment, the first radiation image G1 corresponds to the low-energy image, and the second radiation image G2 corresponds to the high-energy image. Therefore, the attenuation coefficient setting unit 85 derives the attenuation coefficient μls(x,y) of the soft part of the low-energy image and the attenuation coefficient μhs(x,y) of the soft part of the high-energy image by Expression (14) and Expression (15). Note that in Expression (14) and Expression (15), the description of (x,y) is omitted. In addition, μlm is the attenuation coefficient of the muscle in the low-energy image, μlf is the attenuation coefficient of the fat in the low-energy image, μhm is the attenuation coefficient of the muscle in the high-energy image, and μhf is the attenuation coefficient of the fat in the high-energy image.

μls=(1−R)×μlm+R×μlf  (14)

μhs=(1−R)×μhm+R×μhf  (15)

In the fifth embodiment, the subtraction unit 63 performs the subtraction processing by using the attenuation coefficients μls and μhs set by the attenuation coefficient setting unit 85. In this case, the subtraction unit 63 derives the weighting coefficients α and β in Expression (1) and Expression (2). The weighting coefficient α used in the fifth embodiment is derived by α=μls/μhs by using the attenuation coefficient set by the attenuation coefficient setting unit 85. Note that in the fifth embodiment, the attenuation coefficients μlb and μhb of the bone part are used, which are prepared in advance and stored in the storage 53. Further, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by Expression (1) and Expression (2) by using the derived weighting coefficients α and β.

Then, an information derivation device according to a sixth embodiment of the present disclosure will be described. Note that a configuration of the information derivation device according to the sixth embodiment is the same as the configuration of the information derivation device 50 according to the first embodiment, and only the processing performed by the scattered ray removal unit 62 is different, so that the detailed description will be omitted here. FIG. 28 is a diagram showing a functional configuration of the scattered ray removal unit of the information derivation device according to the sixth embodiment. As shown in FIG. 28, a scattered ray removal unit 62A of the information derivation device according to the sixth embodiment comprises an imaging condition acquisition unit 91, a body thickness derivation unit 92, a characteristic acquisition unit 93, a ray distribution derivation unit 94, and a calculation unit 95.

Note that in the sixth embodiment, the first and second radiation images G1 and G2 are acquired by imaging the subject H lying on an imaging table 4 by using an imaging apparatus 1A shown in FIG. 29 by further using a scattered ray removal grid 8, and the first and second radiation images G1 and G2 are used as the teacher data 40. Note that in FIG. 29, the imaging table 4 has a top plate 4A, and the scattered ray removal grid (hereinafter, simply referred to as a grid) 8, the first radiation detector 5, the radiation energy conversion filter 7, and the second radiation detector 6 are disposed below the top plate 4A in order from the side of the radiation source 3, and these components are attached below the top plate 4A by an attachment portion 4B.

In a case of generating the teacher data 40, the imaging condition acquisition unit 91 acquires the imaging conditions used at the time of the energy subtraction imaging of the subject H from the image storage system 9 in order to acquire the first and second radiation images G1 and G2 as the learning data 41. Note that in the sixth embodiment, a radiation quality of the radiation is also included as the imaging condition. The radiation quality is defined by using one or more of the tube voltage [kV] of a radiation generator in the radiation source 3, a total filtration amount [mmAl equivalent], or a half-value layer [mmAl]. The tube voltage means the maximum value of the generated radiation energy distribution. The total filtration amount is obtained by converting the filtration amount of each constituting component which configures the imaging apparatus 1, such as a radiation generator and a collimator, in the radiation source 3 into a thickness of the aluminum. The influence of the beam hardening in the imaging apparatus 1 is larger and the total amount of high-energy components in the wavelength distribution of the radiation is larger as the total filtration amount is larger. The half-value layer is defined by the thickness of the aluminum necessary to attenuate the dose in half with respect to the generated radiation energy distribution. The high-energy components in the wavelength distribution of the radiation is larger as the aluminum in the half-value layer is thicker.

The body thickness derivation unit 92 derives the body thickness distribution of the subject H based on at least one of the first radiation image G1 or the second radiation image G2 and the imaging conditions. Hereinafter, the body thickness distribution derived by the body thickness derivation unit 92 is referred to as the initial body thickness distribution T0. Hereinafter, the derivation of the initial body thickness distribution T0 will be described. Note that in the following description, the derivation of the body thickness and the removal of the scattered ray by using the first radiation image G1 will be described, but the derivation of the body thickness and the removal of the scattered ray can be similarly performed by using the second radiation image G2.

First, in a case in which the radiation source 3 is driven to emit the radiation to the radiation detector 5 in a state in which the subject H is not present, a reaching dose I0(x,y) of the radiation emitted from the radiation source 3 to the radiation detector 5 is represented by Expression (16). In Expression (16), mAs included in the imaging conditions is a tube current-time product, and kV is the tube voltage. Note that the half-value layer is also taken into consideration, the reaching dose I0(x,y) is represented by Expression (16-1). Here, F is a non-linear function that represents the radiation dose that reaches to the radiation detector 5 in a case in which the dose (for example, 1mAs), which is a standard, is emitted to the radiation detector 5 at the SID (for example, 100 cm), which is a standard, in a state in which the subject H is not present. F is changed for each tube voltage or depending on the tube voltage and the half-value layer. In addition, since the reaching dose I0 is derived for each pixel of the radiation image acquired by the radiation detector 5, (x,y) represents the pixel position of each pixel. In addition, in the following description, in order to include both a case in which the half-value layer is considered and a case in which the half-value layer is not considered, each expression is represented by including mmAl in parentheses as shown in Expression (16-2).

I0(x,y)=mAs×F(kV)/SID ²  (16)

I0(x,y)=mAs×F(kV,mmAl)/SID ²  (16-1)

I0(x,y)=mAs×F(kV(mmAl))/SID ²  (16-2)

In addition, in a case in which the initial body thickness distribution is defined as T0, the attenuation coefficient of the subject H in a case of having the initial body thickness distribution T0 is defined as μ(T0), and a scattered-to-primary ratio, which is a ratio between the scattered ray dose included in the radiation after being transmitted through the subject H having the initial body thickness distribution T0 and a primary ray dose in a case in which a scattered ray spread is not considered is defined as STPR(T0), the dose I1 after being transmitted through the subject H is represented by Expression (17). Note that in Expression (17), the initial body thickness distribution T0, the reaching dose I0, and the dose I1 are derived for each pixel of the simple radiation image G0, but (x,y) is omitted. In addition, STPR is a non-linear function that depends on the tube voltage [kV] and the half-value layer [mmAl] in addition to the body thickness, but in Expression (17), kV and mmAl are omitted.

I1=I0×exp{−μ(T0)×T0}×{1+STPR(T0)}  (17)

In Expression (17), the dose I1 is a pixel value in each pixel of the simple radiation image G0, and the reaching dose I0 is derived by Expression (16) and Expression (16-1). On the other hand, since F and STPR are non-linear functions, Expression (17) cannot be algebraically solved for T0. Therefore, the body thickness derivation unit 92 defines an error function E1 shown in Expression (18) or Expression (18-1). Further, T0 at which the error function E1 is minimized or the error function E1 is smaller than a predetermined threshold value is derived as the initial body thickness distribution. In this case, the body thickness derivation unit 92 derives the initial body thickness distribution T0 by using an optimization algorithm, such as the steepest descent method and the conjugate gradient method.

E1=[I1−I0×exp{−μ(T0)×T0}×{1+STPR(T0)}]²  (18)

E1=|I1−I0×exp{−μ(T0)×T0}×{1+STPR(T0)}|  (18-1)

The characteristic acquisition unit 93 acquires a radiation characteristic of an object interposed between the subject H and the first and second radiation detectors 5 and 6 at the time of imaging. Here, in a case in which the radiation after being transmitted through the subject H is transmitted through the object interposed between the subject H and the radiation detector 5, a transmittance of the radiation is changed depending on the radiation quality of the radiation after being transmitted through the subject H. In addition, a primary ray transmittance and a scattered ray transmittance included in the radiation after being transmitted through the subject H are different due to the difference in the traveling direction of the radiation and the quality of the radiation. Therefore, in the sixth embodiment, as the radiation characteristics of the object, the primary ray transmittance and the scattered ray transmittance of the object are used.

Note that, as described above, in a case in which the radiation after being transmitted through the subject H is transmitted through the object interposed between the subject H and the radiation detector 5, a transmittance of the radiation is changed depending on the radiation quality of the radiation after being transmitted through the subject H. In addition, the radiation quality of the radiation after being transmitted through the subject H depends on the body thickness distribution T of the subject H. Therefore, the primary ray transmittance and the scattered ray transmittance can be represented by Tp(T) and Ts(T), respectively, as functions of the body thickness distribution T of the subject H.

In addition, the radiation quality of the radiation after being transmitted through the subject H also depends on the radiation quality of the radiation source 3 included in the imaging conditions. The radiation quality depends on the tube voltage and the half-value layer. Therefore, strictly speaking, the primary ray transmittance and the scattered ray transmittance are represented by Tp(kV(mmAl), T) and Ts(kV(mmAl), T), respectively. Note that in the following description, the primary ray transmittance and the scattered ray transmittance may be simply represented by Tp and Ts.

Here, as described above, the primary ray transmittance Tp and the scattered ray transmittance Ts of the object interposed between the subject H and the radiation detector 5 depend on the body thickness distribution T of the subject H. Therefore, in the sixth embodiment, the primary ray transmittance Tp and the scattered ray transmittance Ts of the object depending on the body thickness distribution T of the subject H are measured by using phantom with various thicknesses that imitate the body thickness distribution T of the subject, a table that defines a relationship between the body thickness distribution T of the subject H and the primary ray transmittance Tp and the scattered ray transmittance Ts of the object is generated based on a result of measurement and stored in the storage of the information derivation device according to the sixth embodiment. Hereinafter, the measurement of the primary ray transmittance Tp and the scattered ray transmittance Ts of the object depending on the body thickness distribution T of the subject H will be described.

First, the calculation of the scattered ray transmittance Ts will be described. FIGS. 30 and 31 are diagrams for describing the measurement of the scattered ray transmittance Ts depending on the body thickness of the subject H. First, as shown in FIG. 30, a phantom 101 that imitates the human body is placed on a surface of the radiation detector 5, and a lead plate 102 is further placed on the phantom 101. Here, the phantom 101 has various thicknesses, such as 5 cm, 10 cm, and 20 cm, and is made of a material, such as acrylic, having a radiation transmittance similar to that of water, for example. In this state, by driving the radiation source 3 and irradiating the radiation detector 5 with the radiation, the characteristic acquisition unit 93 acquires a radiation image K0 for measurement. A signal value of the radiation image K0 is larger in the region in which the radiation is directly emitted to the radiation detector 5, and the signal value is smaller in the order of the region of the phantom 101 and the region of the lead plate 102.

Note that since the lead plate 102 does not transmit the radiation, the signal value should be 0 in a region of the lead plate 102 in the radiation image K0. However, the radiation scattered by the phantom 101 reaches a region corresponding to the lead plate 102 of the radiation detector 5. Therefore, the region of the lead plate 102 in the radiation image K0 has a signal value S0 corresponding to the scattered ray component by the phantom 101.

Next, as shown in FIG. 31, the phantom 101 is placed on the top plate 4A of the imaging apparatus 1, and the lead plate 102 is further placed on the phantom 101. Further, as in a case of imaging the subject H, the characteristic acquisition unit 93 acquires a radiation image K1 for measurement by driving the radiation source 3 to irradiate the radiation detector 5 with the radiation in a state in which the radiation detector 5 and the grid 8 are disposed below the top plate 4A. Similar to the radiation image K0, a signal value of the radiation image K1 is larger in the region in which the radiation is directly emitted to the radiation detector 5, and the signal value is smaller in the order of the region of the phantom 101 and the region of the lead plate 102. Here, as shown in FIG. 31, in a case in which imaging is performed in a state in which the top plate 4A and the grid 8 are interposed between the phantom 101 and the radiation detector 5, the radiation scattered by the top plate 4A and the grid 8 also reaches the region corresponding to the lead plate 102 of the radiation detector 5 in addition to the radiation scattered by the phantom 101. Therefore, the region of the lead plate 102 in the radiation image K1 has a signal value S1 corresponding to the scattered ray component by the phantom 101, the top plate 4A, and the grid 8.

Note that since the signal value S1 includes the scattered ray component due to the top plate 4A and the grid 8, the signal value S1 is larger than the signal value S0 shown in FIG. 30. Therefore, in a case of imaging the phantom 101 having a thickness of t, the scattered ray transmittance Ts of the object interposed between the subject H and the radiation detector 5, that is, the top plate 4A and the grid 8 can be calculated by S1/S0.

In the sixth embodiment, the characteristic acquisition unit 93 calculates the scattered ray transmittance Ts corresponding to each thickness as shown in FIGS. 30 and 31 by using at least two types of phantoms having different thicknesses. In addition, the characteristic acquisition unit 93 derives the scattered ray transmittance Ts having a thickness that is not present in the phantom 101 by interpolating the scattered ray transmittance Ts for a plurality of measured thicknesses. As a result, as shown in FIG. 32, the characteristic acquisition unit 93 generates a table LUT8 that represents a relationship between the body thickness of the subject H and the scattered ray transmittance Ts of the object interposed between the subject H and the radiation detector 5 by interpolating the scattered ray transmittance for the thickness between the thicknesses.

Next, the calculation of the primary ray transmittance will be described. FIGS. 33 and 34 are diagrams for describing the measurement of the primary ray transmittance Tp depending on the body thickness of the subject H. First, as shown in FIG. 33, the phantom 101 that imitates the human body is first placed on the surface of the radiation detector 5. Here, as the phantom 101, the same phantom as in a case in which the scattered ray transmittance Ts is derived is used. Further, in this state, by driving the radiation source 3 and irradiating the radiation detector 5 with the radiation, the characteristic acquisition unit 93 acquires a radiation image K2 for measurement. A signal value S2 in a region corresponding to the phantom 101 in the radiation image K2 includes both the primary ray component and the scattered ray component of the radiation transmitted through the phantom 101. Here, the scattered ray component of the radiation transmitted through the phantom 101 is the signal value S0 in the radiation image K0 obtained by the method shown in FIG. 30. Therefore, the primary ray component of the radiation transmitted through the phantom 101 is derived by S2−S0.

Next, as shown in FIG. 34, the phantom 101 is placed on the top plate 4A of the imaging apparatus 1, and as in a case of imaging the subject H, the characteristic acquisition unit 93 acquires a radiation image K3 for measurement by driving the radiation source 3 to irradiate the radiation detector 5 with the radiation in a state in which the radiation detector 5 and the grid 8 are disposed below the top plate 4A. A signal value S3 in a region corresponding to the phantom 101 in the radiation image K3 includes both the primary ray component and the scattered ray component of the radiation transmitted through the phantom 101, the top plate 4A, and the grid 8. Here, the scattered ray component of the radiation transmitted through the phantom 101, the top plate 4A, and the grid 8 is the signal value S1 in the radiation image K1 obtained by the method shown in FIG. 31. Therefore, the primary ray component of the radiation transmitted through the phantom 101, the top plate 4A, and the grid 8 is derived by S3−S1.

Therefore, it is possible to calculate, by (S3−S1)/(S2−S0), the primary ray transmittance Tp of the top plate 4A and the grid 8 interposed between the subject H and the radiation detector 5 in a case in which the phantom 101 is imaged. Further, in the sixth embodiment, the characteristic acquisition unit 93 calculates the primary ray transmittance Tp corresponding to each thickness as shown in FIGS. 33 and 34 by using at least two types of phantoms having different thicknesses. In addition, the characteristic acquisition unit 93 derives the primary ray transmittance Tp having a thickness that is not present in the phantom 101 by interpolating the primary ray transmittance Tp for a plurality of measured thicknesses. As a result, as shown in FIG. 35, the characteristic acquisition unit 93 generates a table LUT9 that represents a relationship between the body thickness of the subject H and the primary ray transmittance Tp of the object interposed between the subject H and the radiation detector 5.

The tables LUT8 and LUT9, which are generated as described above, are stored in the storage 53 of the information derivation device according to the sixth embodiment. Note that the table is generated depending on various imaging conditions (that is, the radiation quality, the dose, and the radiation source distance), and the type of grid 8 to be used, and is stored in the storage 13.

The characteristic acquisition unit 93 acquires the primary ray transmittance Tp(T0) and the scattered ray transmittance Ts(T0) corresponding to the initial body thickness distribution T0 for the object interposed between the subject H and the radiation detector 5 with reference to the tables LUT8 and LUT9 stored in the storage 53 depending on the imaging conditions acquired by the imaging condition acquisition unit 91. Note that since the primary ray transmittance Tp and the scattered ray transmittance Ts also depend on the radiation quality, the primary ray transmittance Tp and the scattered ray transmittance Ts can be represented by Tp(kV(mmAl), T0) and Ts(kV(mmAl), T0), respectively.

The ray distribution derivation unit 94 derives the primary ray distribution and the scattered ray distribution of the radiation detected by the radiation detector 5 by using the imaging conditions, the body thickness distribution, and the radiation characteristic of the object interposed between the subject H and the radiation detector 5. Here, a primary ray distribution Ip0 and a scattered ray distribution Is0 after being transmitted through the subject H are represented by Expression (19) and Expression (20) in a case in which the body thickness distribution is defined as T. PSF in Expression (20) is a point spread function that represents the distribution of the scattered rays spreading from one pixel, and is defined depending on the radiation quality and the body thickness. In addition, * indicates a convolution operation. The primary ray distribution Ip0 and the scattered ray distribution Is0 are derived for each pixel of the simple radiation image G0, but (x,y) are omitted in Expression (19) and Expression (20). In addition, in the sixth embodiment, derivation of the body thickness distribution, the primary ray distribution Ip0, and the scattered ray distribution Is0 is repeatedly performed as described below, but the first primary ray distribution Ip0 and the scattered ray distribution Is0 At the time of the derivation, the initial body thickness distribution T0 is used as the body thickness distribution T.

Ip0=I0×exp{−μ(T)×T}  (19)

Is0=Ip0×STPR(kV(mmAl),T)*PSF(kV(mmAl),T)  (20)

Moreover, the ray distribution derivation unit 94 derives a primary ray distribution Ip1 and a scattered ray distribution Is1 reaching to the radiation detector 5 by Expression (21) and Expression (22) by using the primary ray transmittance Tp and the scattered ray transmittance Ts of the object interposed between the subject H and the radiation detector 5. Moreover, a sum Iw1 of the primary ray distribution Ip1 and the scattered ray distribution Is1 is derived by Expression (23). Also in Expression (21) and Expression (22), the initial body thickness distribution T0 is used as the body thickness distribution T in a case of first derivation of the first primary ray distribution Ip1 and the scattered ray distribution Is1.

Ip1=Ip0×Tp(kV(mmAl),T)  (21)

Is1=Is0×Ts(kV(mmAl),T)  (22)

Iw1=Ip1+Is1  (23)

The calculation unit 95 derives an error E2 between the sum Iw1 of the primary ray distribution Ip1 and the scattered ray distribution Is1 and the dose at each pixel position of the first radiation image G1, that is, a pixel value I1. The derivation of the error E2 is performed by Expression (24) or Expression (24-1). In Expression (24) and Expression (24-1), N represents the number of pixels of the first radiation image G1 and Σ represents a sum of all of the pixels of the first radiation image G1. Note that since in Expression (24-1), I1/Iw1 is calculated in the log, the error E2 can be derived without depending on the emitted dose to the subject H, that is, the reaching dose 10.

E2=(1/N)×Σ{I1−Iw1}²  (24)

E2=(1/N)×Σ|log{I1/Iw1}|  (24-1)

Further, the calculation unit 95 updates the body thickness distribution T such that the error E2 is minimized or the error E2 is smaller than a predetermined threshold value. Further, the calculation unit 95 repeats the acquisition of the primary ray transmittance Tp and the scattered ray transmittance Ts based on the updated body thickness distribution, and the derivation of the primary ray distribution Ip1 and the scattered ray distribution Is1. Here, the calculation performed by the calculation unit 95 is referred to as a repetition calculation. In addition, in the sixth embodiment, the calculation unit 95 repeatedly performs the repetition calculation such that the error E2 is smaller than the predetermined threshold value. Further, the calculation unit 95 outputs a processed first radiation image Gc1 which uses, as the pixel value, the primary ray distribution Ipc derived based on the body thickness distribution Tc of the subject H in which the error E2 is smaller than the predetermined threshold value.

Note that the characteristic acquisition unit 93 and the ray distribution derivation unit 94 perform the repetition acquisition of the primary ray transmittance Tp and the scattered ray transmittance Ts and the repetition derivation of the primary ray distribution Ip1 and the scattered ray distribution Is1, respectively.

On the other hand, for the second radiation image G2, the primary ray distribution Ipc is derived in the same manner as for the first radiation image G1. Note that the primary ray distribution for the first radiation image G1 is defined as Ipc-1, and the primary ray distribution for the second radiation image G2 is defined as Ipc-2. Then, the calculation unit 95 outputs a processed second radiation image Gc2 which has the primary ray distribution Ipc-2 as the pixel value.

In the sixth embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by using the processed first and second radiation images Gc1 and Gc2.

Note that in a case in which the simple radiation image G0 obtained by simply imaging the subject H by the imaging apparatus 1A shown in FIG. 29 is input, the learned neural network 23A constructed by performing the learning by using the bone part image Gb and the soft part image Gs derived in the sixth embodiment as the correct answer data outputs the bone part image Gb and the soft part image Gs in which the object interposed between the subject H and the radiation detector is considered.

In addition, in the sixth embodiment, the top plate 4A and the grid 8 of the imaging table 4 are used as the objects interposed between the subject H and the radiation detector 5, but as shown in FIG. 36, an air layer 103 may be interposed between the top plate 4A and the grid 8. In such a case, it is preferable that the ray distribution derivation unit 94 derive the primary ray distribution Ip1 and the scattered ray distribution Is1 by including the air layer 103 as the object interposed between the subject H and the radiation detector 5. In this case, as shown in Expression (21-1) and Expression (22-1), the primary ray distribution Ip1 and the scattered ray distribution Is1 need only be derived by performing convolution operation on a point spread function PSFair(kV(mmAl), tair) depending on the thickness tair of the air layer 103 with respect to Expression (21) and Expression (22). Note that the thickness tair of the air layer 103 is the distance between the lower surface of the top plate 4A and the surface of the grid 8 on the subject H side.

Ip1=Ip0×Tp(kV(mmAl),T)*PSF_(air)(kV(mmAl),t _(air))  (21-1)

Is1=Is0×Ts(kV(mmAl),T)*PSF_(air)(kV(mmAl),t _(air))  (22-1)

Then, an information derivation device according to a seventh embodiment of the present disclosure will be described. Note that a configuration of the information derivation device according to the seventh embodiment is the same as the configuration of the information derivation device according to the first embodiment, and only the processing performed by the scattered ray removal unit 62 is different, so that the detailed description will be omitted here. FIG. 37 is a diagram showing a functional configuration of the scattered ray removal unit of the information derivation device according to the seventh embodiment. As shown in FIG. 37, a scattered ray removal unit 62B of the information derivation device according to the seventh embodiment comprises a first derivation unit 97, a second derivation unit 98, and an image generation unit 99.

The information derivation device according to the seventh embodiment performs the scattered ray removal processing by causing the first derivation unit 97 to derive the first primary ray distribution and the scattered ray distribution of the radiation transmitted through the subject H by using the first and second radiation images G1 and G2, causing the second derivation unit 98 to derive the second primary ray distribution and the scattered ray distribution of the radiation transmitted through the object by using the first primary ray distribution, the scattered ray distribution, and the radiation characteristic of the object interposed between the subject H and the radiation detector that detects the radiation image, and causing the image generation unit 99 to derive the radiation images after transmission through the subject and the object by using the second primary ray distribution and the scattered ray distribution. Note that the scattered ray removal processing performed by the information derivation device according to the seventh embodiment is disclosed in WO2020/241664A. Hereinafter, the scattered ray removal processing from the first radiation image G1 will be described, but the scattered ray removal processing can also be performed on the second radiation image G2 in the same manner.

The first derivation unit 97 uses the first radiation image G1 to estimate the components (the primary ray component and the scattered ray component) of the radiation transmitted through the subject H. In derivation processing performed by the first derivation unit 97 (hereinafter referred to as first derivation processing), the “radiation transmitted through the subject H” refers to the radiation transmitted through the subject H after being transmitted through the subject H and before transmitted through the object, such as the top plate 4A and the grid 8.

In addition, the component of the radiation transmitted through the subject H specifically means the component of the radiation transmitted through the subject H and/or the component of the radiation scattered by the subject H. That is, the component of the radiation transmitted through the subject H is the primary ray component after being transmitted through the subject H. The component of the radiation scattered by the subject H is the scattered ray component after being transmitted through the subject H. Regarding the radiation incident on the subject H toward any position X, in a case in which the subject H is regarded as an operator g1 that generates the primary ray component and an operator h1 that generates the scattered ray component, as shown in FIG. 38, the primary ray component after being transmitted through the subject H is g1(X) and the scattered ray component after being transmitted through the subject H is h1(X).

The first derivation unit 97 uses the first radiation image G1 to estimate the primary ray component g1(X) after being transmitted through the subject H, the scattered ray component h1(X) after being transmitted through the subject H, or both the primary ray component g1(X) and the scattered ray component h1(X). In the seventh embodiment, the first derivation unit 97 uses the first radiation image G1 to estimate each of the primary ray component g1(X) after being transmitted through the subject H and the scattered ray component h1(X) after being transmitted through the subject H.

Note that in a case in which the primary ray component g1(X) after being transmitted through the subject H is estimated from the first radiation image G1, the first derivation unit 97 estimates the scattered ray component h1(X) after transmitted through the subject H by subtracting the estimated primary ray component g1(X) from the first radiation image G1. In addition, in a case in which the scattered ray component h1(X) after being transmitted through the subject H is estimated from the first radiation image G1, the first derivation unit 97 estimates the primary ray component g1(X) after transmitted through the subject H by subtracting the estimated scattered ray component h1(X) from the first radiation image G1.

The first derivation processing performed by the first derivation unit 97 can be performed by, for example, estimating the body thickness of the subject H by using the first radiation image G1, and estimating the components of the radiation transmitted through the subject H by using the estimated body thickness of the subject H. In this case, the first derivation unit 97 estimates the primary ray component g1(X) of the radiation transmitted through the subject H and the scattered ray component h1(X) of the radiation scattered by the subject H for each pixel of the first radiation image G1 (or for each predetermined compartment consisting of a plurality of pixels) based on the estimated body thickness of the subject H.

For example, as shown in FIG. 39, a pixel value V2 in a case in which the subject H is present (“with the subject”) is smaller than a pixel value V1 in a case in which the subject H is not present (“without the subject”). This is due to absorption of the subject H and the like. Therefore, a difference Δ(=V1−V2) between the pixels values relates to the body thickness of the subject H. On the other hand, the pixel value V1 in a case in which the subject H is not present can be known by, such as the pixel value of the region (direct region) in which the radiation reaches the radiation detector 5 without being transmitted through the subject H, or the experiment performed in advance (imaging without the subject H). Therefore, the first derivation unit 97 can estimate the body thickness of the subject H from the pixel value V2 of the first radiation image G1 imaged with the subject H.

In addition, the primary ray component g1(X) and the scattered ray component h1(X) after being transmitted through the subject H both relate to the body thickness of the subject H. For example, as the body thickness of the subject H is thicker, the primary ray component g1(X) is smaller due to the absorption of the subject H and the like, and the scattered ray component h1(X) is larger with respect to the incident radiation. Such properties of the subject H, that is, a transmission amount and a scattering amount of the subject H with respect to the radiation having a specific energy can be obtained in advance by an experiment or the like before radiography.

For this reason, the first derivation unit 97 holds the characteristics relating to the transmission amount and the scattering amount (hereinafter, referred to as the subject scattering characteristics) in the form of a function or a table or the like, for example, for each subject H or for each imaging part of the subject H. Further, by using the energy of the radiation used for imaging and the estimated actual body thickness of the subject H to obtain the transmission amount and the scattering amount of the radiation, the primary ray component g1(X) and the scattered ray component h1(X) after being transmitted through the subject H are estimated.

Results of derivation output by the first derivation unit 97 (hereinafter referred to as the first result of derivation) are the primary ray component g1(X) after being transmitted through the subject H at a position P1, the scattered ray component h1(X) after being transmitted through the subject H at the position P1, or an intensity distribution f1(X) of the radiation after being transmitted through the subject H at the position P1. The intensity distribution f1(X) of the radiation at the position P1 is, for example, the sum or weighted sum of the primary ray component g1(X) and the scattered ray component h1(X). In the seventh embodiment, the first derivation unit 97 outputs the intensity distribution f1(X) of the radiation after being transmitted through the subject H at the position P1 as a first result of derivation in the form of an image or the form of a collection of data capable of constructing an image, for example. Note that the first derivation unit 97 can also output one of the primary ray component g1(X) or the scattered ray component h1(X) after being transmitted through the subject H as the result of derivation.

The second derivation unit 98 estimates the components of the radiation transmitted through the object by using the result of derivation of the first derivation unit 97 and the scattering characteristics of the object through which the radiation transmitted through the subject H is further transmitted, such as the top plate 4A and the grid 8. In derivation processing of the second derivation unit 98 (hereinafter referred to as second derivation processing), “transmitted through the object” means that the radiation is transmitted through a certain position in which the subject H is present and then transmitted through the object. Therefore, depending on a specific shape of the subject H, the radiation is directly transmitted through the object without being transmitted through the subject H.

Specifically, the second derivation unit 98 estimates the component of the radiation transmitted through the subject H and the object, or the component of the radiation scattered by at least one of the subject H or the object. The component of the radiation transmitted through the subject H and the object is the primary ray component after being transmitted through the object. The radiation component scattered by at least one of the subject H or the object is the scattered ray component after being transmitted through the object.

The scattering characteristics of the object determine the distribution of the radiation doses transmitted through the object and/or the radiation dose scattered by the object. In the seventh embodiment, a scattering characteristic f2(X) including a first characteristic g2(X) that determines the distribution of the radiation doses transmitted through the object and a second characteristic h2(X) that determines the distribution of the radiation doses scattered by the object is used. Specifically, the scattering characteristic f2(X) is the sum or weighted sum of the first characteristic g2(X) and the second characteristic h2(X), for example, f2(X)=g2(X)+h2(X).

The first characteristic g2(X) is a function or a table or the like that determines the transmitted dose of the radiation directly incident on the object toward any position X without being transmitted through the subject H. In addition, the second characteristic h2(X) is a function or a table or the like that determines the transmitted dose of the radiation directly incident on the object toward any position X without being transmitted through the subject H. For example, in a case in which the object is only the top plate 4A of the imaging table 4, the first characteristic g2(X) determines the distribution of the transmitted doses of the top plate 4A, and the second characteristic h2(X) determines the distribution of the scattered ray doses of the top plate 4A. A state of a specific configuration of the object (use or non-use of the imaging table 4 or the like) is known before radiography. Therefore, the first characteristic g2(X) and the second characteristic h2(X) can be obtained in advance by experiments or the like, for example, for each specific configuration of the object or for each combination of the objects. In addition, in a case in which the object is regarded to generate the primary ray component and the scattered ray component from the incident radiation, the first characteristic g2(X) is the operator that generates the primary ray component depending on the incident radiation, and the second characteristic h2(X) is the operator that generates the scattered ray component depending on the incident radiation.

In the seventh embodiment, the second derivation unit 98 holds the first characteristic g2(X) and the second characteristic h2(X) in advance for each specific configuration of the object, for example. As a result, the second derivation unit 98 holds the scattering characteristic f2(X) of the object in advance. Note that the second derivation unit 98 can acquire the first characteristic g2(X), the second characteristic h2(X), and/or the scattering characteristic f2(X), as needed.

As shown in FIG. 40, the intensity distribution after the radiation (X) incident toward any point X0 is transmitted through the object can be approximated by a point spread function (PSF) 120. For example, the PSF 120 is a Gaussian function. Further, of the radiation (X) incident on the object toward any point X0, the component reaching any point X0 and its vicinity is a distribution 121 of the primary ray component, and a portion excluding the distribution 121 of the primary ray component from the PSF 120 is a distribution 122 of the scattered ray component.

Moreover, since the energy of the radiation used for imaging, the material (density and the like), and the thickness (mass) of the top plate 4A, which is the object, are known, the specific shape of the PSF 120, such as peak height and full width at half maximum, is predetermined. Therefore, for example, the second characteristic h2(X) can be obtained in advance by performing a deconvolution operation of the distribution 122 of the scattered ray component on the radiation image obtained by performing imaging without the subject H. In addition, the first characteristic g2(X) can be obtained in advance by subtracting the second characteristic h2(X) from the radiation image obtained by performing imaging without the subject H, or by performing the deconvolution operation of the distribution 121 of the primary ray component.

The second derivation unit 98 estimates the component of the radiation transmitted through the object by applying the scattering characteristic of the object to the first result of derivation which is the result of derivation of the first derivation unit 97. Specifically, the object receives the incidence of the radiation of the distribution represented by the first result of derivation. Therefore, the second derivation unit 98 estimates the component of the radiation transmitted through the object by setting an argument of the scattering characteristic f2(X) of the object as the first result of derivation (f1(X)). That is, the second derivation unit 98 estimates the component of the radiation transmitted through the object by the calculation based on Expression (25). In the seventh embodiment, since the first result of derivation is f1(X)=g1(X)+h1(X), Expression (25) can be represented as Expression (26), and can be expanded as in Expression (27).

f2(f1(X))=g2(f1(X))+h2(f1(X))  (25)

f2(f1(X))=g2(g1(X)+h1(X))+h2(g1(X)+h1(X))  (26)

f2(f1(X))=g2g1(X)+g2h1(X)+h2g1(X)+h2h1(X)  (27)

As shown in FIG. 41, the first term “g2g1(X)” on the right side of Expression (27) represents radiation Ra1 of the radiation used for imaging, which is transmitted through the subject H and transmitted through the object (top plate 4A in FIG. 41), and then reaches a pixel P(X0) at any point X0. The second term “g2h1(X)” on the right side of Expression (27) represents radiation Ra2 of the radiation used for imaging, which is scattered by a scattering body D1 included in the subject H and then is transmitted through the object to reach the pixel P(X0) at any point X0. The third term “h2g1(X)” on the right side of Expression (27) represents radiation Ra3 of the radiation used for imaging, which is transmitted through the subject H and then is scattered by a scattering body D3 included in the object to reach the pixel P(X0) at any point X0. In addition, the fourth term “h2h1(X)” on the right side of Expression (27) represents radiation Ra4 of the radiation used for imaging, which is scattered by a scattering body D2 included in the subject H and then is further scattered by a scattering body D4 included in the object to reach the pixel P(X0) at any point X0.

As described above, the second derivation unit 98 obtains the first term “g2g1(X)” of Expression (27) and/or the sum of the second to fourth terms “g2h1(X)+h2g1(X)+h2h1(X)”.

This is because the first term “g2g1(X)” of Expression (27) represents the distribution of the primary ray components after being transmitted through the object, and the sum of the second to fourth terms “g2h1(X)+h2g1(X)+h2h1(X)” represents the distribution of the scattered ray components after being transmitted through the object. In the seventh embodiment, the second derivation unit 98 obtains the distribution g2g1(X) of the primary ray components after being transmitted through the object, and outputs the obtained distribution g2g1(X) as the result of derivation.

The image generation unit 99 uses the result of derivation of the second derivation unit 98 to generate the processed radiation image that forms the image of the subject H by the radiation transmitted through the subject H and the object. In a case in which the second derivation unit 98 estimates the primary ray component of the radiation transmitted through the subject H and the object, the image generation unit 99 images a second result of derivation, which is the result of derivation of the second derivation unit 98, to generate the processed radiation image. In addition, in a case in which the second derivation unit 98 estimates the scattered ray component of the radiation scattered by the subject H or the object, the image generation unit 99 subtracts the second result of derivation, which is the result of derivation of the second derivation unit 98, from the first radiation image G1 to generate the processed radiation image.

In the seventh embodiment, since the second derivation unit 98 outputs the distribution of the primary ray component after being transmitted through the object, the image generation unit 99 generates the processed radiation image by imaging the output distribution of the primary ray component. Therefore, the distribution g2g1(X) of the primary ray component output by the second derivation unit 98 is the processed radiation image in which the scattered ray component is substantially removed. Note that the image generation unit 99 can perform various pieces of image processing (for example, contrast adjustment processing or structure emphasis processing) on the generated processed radiation image, as needed.

Note that the scattered ray removal processing according to the sixth and seventh embodiments may be performed in the first to fifth embodiments. In this case, the information derivation devices 50A to 50D include the scattered ray removal unit 62A or the scattered ray removal unit 62B instead of the scattered ray removal unit 62.

Then, an eighth embodiment of the present disclosure will be described. FIG. 42 is a diagram showing a functional configuration of an information derivation device according to the eighth embodiment. Note that in FIG. 42, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. In the seventh embodiment of the present disclosure, the first and second radiation images G1 and G2 are subjected to granulation suppression processing, and the energy subtraction processing is performed by using the first and second radiation images G1 and G2 to which the granulation suppression processing is performed. Therefore, as shown in FIG. 42, an information derivation device 50E according to the eighth embodiment further comprises a processing content derivation unit 86 and a granulation suppression processing unit 87 with respect to the information derivation device 50 according to the first embodiment.

The processing content derivation unit 86 derives a processing content of first granulation suppression processing on the first radiation image G1, and derives a processing content of second granulation suppression processing on the second radiation image G2 based on the processing content of the first granulation suppression processing.

Here, in the present embodiment, the subject H is imaged by the one-shot method to acquire the first and second radiation images G1 and G2. In the case of the one-shot method, the radiation transmitted through the subject H is emitted to the two radiation detectors 5 and 6 stacked with the radiation energy conversion filter 7 interposed therebetween. Therefore, the emitted dose of the second radiation detector 6 on the side away from the radiation source 3 is smaller than that of the first radiation detector 5 on the side closer to the radiation source 3. As a result, the second radiation image G2 has more radiation quantum noise and lower S/N than the first radiation image G1. Therefore, particularly for the second radiation image G2, it is necessary to perform the granulation suppression processing of suppressing granulation caused by quantum noise.

The processing content derivation unit 86 derives the processing content of the first granulation suppression processing on the first radiation image G1 having higher S/N among the first radiation image G1 and the second radiation image G2. Further, the processing content derivation unit 86 derives the processing content of the second granulation suppression processing on the second radiation image G2 based on the processing content of the first granulation suppression processing. Hereinafter, the derivation of the processing content will be described.

Examples of the granulation suppression processing include filtering processing by using a smoothing filter, such as a Gaussian filter, having a predetermined size such as 3×3 or 5×5 centered on an attention pixel. However, in a case in which the Gaussian filter is used, the edges of the structures included in the first and second radiation images G1 and G2 may be blurred. Therefore, in the eighth embodiment, the granulation suppression processing is performed by using an edge storage smoothing filter that suppresses the granulation while preventing edge blurriness. As the edge storage smoothing filter, a bilateral filter, which weights the normal distribution in which the weight of the pixel adjacent to the attention pixel is smaller as the distance from the attention pixel is larger, and the weight is reduced as the difference in the pixel value from the attention pixel is larger, is used.

FIG. 43 is a diagram showing an example of the bilateral filter on the first radiation image G1. Note that FIG. 43 shows two local regions A1 having 5×5 pixels in the vicinity of the edge included in the first radiation image G1 side by side. In addition, the two local regions A1 shown in FIG. 43 are the same, but the positions of the attention pixels are different from each other. In the local region A1 on the left side of FIG. 43, the low-concentration pixel in contact with the boundary of the edge is an attention pixel P11. As described above, the bilateral filter weights the normal distribution in which the weight of the pixel adjacent to the attention pixel is smaller as the distance from the attention pixel is larger, and the weight is reduced as the difference in the pixel value from the attention pixel is larger.

Therefore, the processing content derivation unit 86 determines a filter size of the bilateral filter based on the difference between the pixel value of the attention pixel and the pixel value of the pixels around the attention pixel. For example, as the difference between the pixel value of the attention pixel and the pixel value around the attention pixel is smaller, the filter size is larger. In addition, the processing content derivation unit 86 determines the weight of the bilateral filter based on the difference between the pixel value of the attention pixel and the pixel value of the pixels around the attention pixel. For example, as the difference between the pixel value of the attention pixel and the pixel value around the attention pixel is smaller, the weight of the pixel closer to the attention pixel is larger than the weight of the pixel away from the attention pixel.

As a result, the processing content derivation unit 86 derives a 3×3 bilateral filter F11 having a weight as shown on the left side of FIG. 43 for the attention pixel P11 in the local region A1 on the left side of FIG. 43 as the processing content of the first granulation suppression processing.

In the local region A1 on the right side of FIG. 43, the high-concentration pixel in contact with the boundary of the edge is an attention pixel P12. Therefore, the processing content derivation unit 86 derives a 3×3 bilateral filter F12 having a weight as shown on the right side of FIG. 43 for the attention pixel P12 in the local region A1 on the right side of FIG. 43 as the processing content of the first granulation suppression processing.

FIG. 44 is a diagram showing a local region A2 of the second radiation image corresponding to the local region A1 of the first radiation image shown in FIG. 43. Note that FIG. 44 shows the local region A2 having 5×5 pixels in the vicinity of the edge included in the second radiation image G2. The local region A2 is a region at the same position as the local region A1 in the first radiation image G1 shown in FIG. 43. The two local regions A2 shown in FIG. 44 are the same, but the positions of the attention pixels are different from each other. In the local region A2 on the left side of FIG. 44, the pixel corresponding to the attention pixel P11 shown in the local region A1 on the left side of FIG. 43 is an attention pixel P21. In the local region A2 on the right side of FIG. 44, the pixel corresponding to the attention pixel P12 shown in the local region A1 on the right side of FIG. 43 is an attention pixel P22.

Here, the dose of the radiation emitted to the second radiation detector 6 from which the second radiation image G2 is acquired is lower than the dose of the radiation emitted to the first radiation detector 5 from which the first radiation image G1 is acquired. Therefore, the second radiation image G2 has more radiation quantum noise than the first radiation image G1 and has poor granulation, so that the edges do not clearly appear. In addition, due to the influence of the quantum noise, the low-concentration pixels may be included in the high-concentration region in the vicinity of the boundary of the edge, or the high-concentration pixels may be included in the low-concentration region. Therefore, from the second radiation image G2, it is not possible to appropriately determine the bilateral filter that suppresses granulation while storing the edges, unlike the first radiation image G1.

In this case, it is conceivable to use the smoothing filter, such as the Gaussian filter. However, in a case in which the smoothing filter is used, it is not possible to achieve both noise suppression and edge storage, and as a result, the edges are buried in the noise, and the edges of the structure included in the second radiation image G2 cannot be restored.

Therefore, in the eighth embodiment, the processing content derivation unit 86 derives the processing content of second granulation suppression processing on the second radiation image based on the processing content of the first granulation suppression processing on the first radiation image G1. That is, the processing content derivation unit 86 derives the processing content of the second granulation suppression processing such that the processing content of the second granulation suppression processing performed on each pixel of the second radiation image G2 is the same as the processing content of the first granulation suppression processing performed on the pixel of the first radiation image G1 corresponding to each pixel of the second radiation image G2. Specifically, the processing content derivation unit 86 derives the bilateral filter having the same size and the same weight as the bilateral filter determined for each pixel of the first radiation image G1 as the processing content of the second granulation suppression processing on the second radiation image G2.

FIG. 45 a diagram showing an example of the bilateral filter for the second radiation image G2. Note that FIG. 45 shows the local region A2 having 5×5 pixels in the vicinity of the edge included in the second radiation image G2, as in FIG. 43. As shown in FIG. 45, the processing content derivation unit 86 derives a bilateral filter F21 having the same size and the same weight as the bilateral filter F11 derived for the attention pixel P11 of the local region A1 of the first radiation image G1 as the processing content of the second granulation suppression processing, for the attention pixel P21 of the local region A2 of the second radiation image G2.

In addition, the processing content derivation unit 86 derives a bilateral filter F22 having the same size and the same weight as the bilateral filter F12 derived for the attention pixel P12 of the local region A1 of the first radiation image G1 as the processing content of the second granulation suppression processing, for the attention pixel P22 of the local region A2 of the second radiation image G2.

Note that the processing content of the second granulation suppression processing is derived, it is necessary to associate the pixel positions of the first radiation image G1 and the second radiation image G2 with each other. Therefore, it is preferable to perform registration of the first radiation image G1 and the second radiation image G2.

The granulation suppression processing unit 87 performs the granulation suppression processing on the first radiation image G1 and the second radiation image G2. That is, the granulation suppression processing is performed on the first radiation image G1 and the second radiation image G2 based on the processing content derived by the processing content derivation unit 86. Specifically, the first radiation image G1 is subjected to the filtering processing by the bilateral filter derived for the first radiation image G1. In addition, the second radiation image G2 is subjected to the filtering processing by the bilateral filter derived based on the first radiation image G1.

In the eighth embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by Expression (1) and Expression (2) by using the first and second radiation images G1 and G2 that are subjected to the granulation suppression processing.

Then, a ninth embodiment of the present disclosure will be described. FIG. 46 is a diagram showing a functional configuration of an information derivation device according to the ninth embodiment. Note that in FIG. 46, the same reference numerals are assigned to the same configurations as those in FIG. 42, and the detailed description thereof will be omitted. An information derivation device 50F according to the ninth embodiment of the present disclosure further comprises a map derivation unit 88 that derives a physical quantity map of the subject H based on at least one of the first radiation image G1 or the second radiation image G2 with respect to the information derivation device 50E according to the eighth embodiment, and is different from the eighth embodiment in that the processing content derivation unit 86 derives the processing content of the second granulation suppression processing on the second radiation image G2 based on the physical quantity map.

The map derivation unit 88 derives the physical quantity map for the subject H. Examples of the physical quantity include the body thickness and the bone density of the subject H. As the body thickness, the body thickness distribution in a case in which the scattered ray removal unit 62 satisfies the termination condition can be used. As the bone density, it is possible to use the bone density derived from the bone part image Gb by a predetermined method.

Here, the contrasts of the structures included in the first and second radiation images G1 and G2, which are the targets of the granulation suppression processing, vary depending on the imaging conditions. Therefore, in a case in which the edge storage smoothing processing is performed by using the bilateral filter, it is necessary to control the intensity of the edge to be stored depending on the imaging conditions. On the other hand, in a body thickness map representing the body thickness of the subject H, the contrast of the structure included in the map is represented by the thickness (mm) that does not depend on the imaging conditions.

Therefore, in the ninth embodiment, the processing content derivation unit 86 derives the processing content of the first granulation suppression processing on the first radiation image G1 based on the physical quantity map. FIG. 47 is a diagram showing an example of the bilateral filter for the physical quantity map. Note that FIG. 47 shows two 5×5 pixel local regions A3 in the vicinity of the edges included in the physical quantity map side by side. The two local regions A3 shown in FIG. 47 are the same, but the positions of the attention pixels are different from each other. In the local region A3 on the left side of FIG. 47, the high-concentration pixel on the edge is an attention pixel P31. Therefore, for the attention pixel P31 in the local region A3 on the left side of FIG. 47, a 3×3 bilateral filter F31 having the weight as shown on the left side of FIG. 47 is derived as the processing content of the first granulation suppression processing.

In the local region A3 on the right side of FIG. 47, the low-concentration pixel in contact with the boundary of the edge is an attention pixel P32. Therefore, for the attention pixel P32 in the local region A3 on the right side of FIG. 47, a 3×3 bilateral filter F32 having the weight as shown on the right side of FIG. 47 is derived as the processing content of the first granulation suppression processing.

Also in the ninth embodiment, the processing content derivation unit 86 derives the bilateral filter having the same bilateral filter determined for each pixel of the first radiation image G1 as the processing content of the second granulation suppression processing on the second radiation image G2. That is, in the local region of the second radiation image G2 corresponding to the local region A3 of the first radiation image G1, for the pixel corresponding to the attention pixel P31 in the local region A3, the bilateral filter having the same size and the same weight as the bilateral filter F31 is derived as the processing content of the second granulation suppression processing. In addition, for the pixel of the second radiation image G2 corresponding to the attention pixel P32 of the local region A3 of the first radiation image G1, the bilateral filter having the same size and the same weight as the bilateral filter F32 is derived as the processing content of the second granulation suppression processing.

Further, in the ninth embodiment, the granulation suppression processing is performed on the first radiation image G1 and the second radiation image G2 based on the processing content derived by the processing content derivation unit 86 by the granulation suppression processing unit 87. The subtraction unit 63 derives the bone part image Gb and the soft part image Gs by Expression (1) and Expression (2) by using the first and second radiation images G1 and G2 that are subjected to the granulation suppression processing.

Note that the granulation suppression processing performed in the eighth and ninth embodiments may be performed in the first to seventh embodiments.

Then, an information derivation device according to a tenth embodiment of the present disclosure will be described. Note that a configuration of the information derivation device according to the tenth embodiment is the same as the configuration of the information derivation device according to the first embodiment, and only the processing performed by the scattered ray removal unit 62 is different, so that the detailed description will be omitted here. Note that in the tenth embodiment, the first and second radiation images G1 and G2 acquired by imaging the subject H by using the imaging apparatus 1A shown in FIG. 29 are used as the teacher data 40 as in the sixth embodiment.

FIG. 48 is a diagram showing a functional configuration of the scattered ray removal unit of the information derivation device according to the tenth embodiment. As shown in FIG. 48, a scattered ray removal unit 62C of the information derivation device according to the tenth embodiment comprises a region detection unit 150, a scattered ray image derivation unit 151, a pixel value calculation unit 152, a pixel value substitution unit 153, a boundary position adjustment unit 154, and a calculation unit 155.

The region detection unit 150 obtains the region detection image by detecting a subject region in which the radiation is transmitted through the subject H and reaches the first and second radiation detectors 5 and 6, and a direct radiation region in which the radiation is transmitted through only the top plate 4A and directly reaches the first and second radiation detectors 5 and 6 without being transmitted through the subject H in the first and second radiation images G1 and G2.

Specifically, as shown in FIG. 49, the region detection unit 150 detects the region in which the pixel value of the first radiation image G1 is smaller than a region threshold value as a subject region 160, and detects the region in which the pixel value of the first radiation image G1 is equal to or larger than the region threshold value as a direct radiation region 161. For example, in a case in which the region detection image is a binarized image in which the subject region 160 is “0” and the direct radiation region 161 is “1”, the distribution shown in FIG. 50 is obtained in a specific line 162 of the first radiation image G1. Note that the region detection unit 150 may use the learned neural network that is learned to detect the subject region and the direct radiation region instead of the region detection by threshold value processing.

Note that the region threshold value is determined for each imaging condition, and it is preferable to use the region threshold value corresponding to the imaging condition in a case of imaging the subject H. In addition, region detection processing may be performed on both the first radiation image G1 and the second radiation image G2, but may be performed on any one of the first radiation image G1 or the second radiation image G2.

The scattered ray image acquisition unit 151 obtains the scattered ray image relating to the scattered ray component based on the region detection image and scattered ray spread information relating to the spread of the scattered ray. The calculation unit 155 obtains the radiation image from which the scattered ray component is removed by subtracting the scattered ray image from the first and second radiation images G1 and G2. Note that as the scattered ray spread information used in the scattered ray image acquisition unit 151, for example, the information shown in FIG. 40 can be used.

Further, the scattered ray image derivation unit 151 obtains the scattered ray image including the scattered ray component by the calculation of “region detection image*PSF (* is a convolution operator)” that convolves the region detection image and the PSF. Note that the scattered ray image derivation unit 151 derives the scattered ray image for each of the first radiation image G1 and the second radiation image G2.

In the scattered ray image, as the pixel value of the scattered ray component in the direct radiation region, the pixel value of the direct radiation region 161 of the first and second radiation images G1 and G2 may be used as it is, but the pixel value of the direct radiation region 161 is often saturated (exceeding the pixel value that imaging sensors of the first and second radiation detectors 5 and 6 can receive). Therefore, as the pixel value in the direct radiation region of the scattered ray image, it is preferable theoretically calculate an unsaturated scattered ray pixel value, which is a unsaturated pixel value in which the pixel value is not saturated and the influence of the scattered ray is taken into consideration, and to substitute the theoretically calculated unsaturated scattered ray pixel value with the pixel value in the direct radiation region of the scattered ray image.

The pixel value calculation unit 152 calculates the unsaturated scattered ray pixel value corresponding to the imaging dose by referring to an unsaturated scattered ray pixel value relationship representing a relationship between the dose of the radiation and the unsaturated scattered ray pixel value, which is the unsaturated pixel value in which the pixel value is not saturated and the influence of the scattered ray is taken into consideration. Further, the pixel value substitution unit 153 substitutes the pixel value in the direct radiation region with the unsaturated scattered ray pixel value. Specifically, as shown in FIG. 51, in the scattered ray image derived by the scattered ray image derivation unit 151, the pixel value of the direct radiation region 161 is substituted by the unsaturated scattered ray pixel value obtained by the pixel value calculation unit 152.

The unsaturated scattered ray pixel value relationship is predetermined in a table (not shown) for the direct radiation region for each imaging condition, and is stored in the storage 53. The pixel value calculation unit 152 refers to the table for the direct radiation region to use any of the unsaturated scattered ray pixel value relationships corresponding the imaging condition in a case of imaging the subject H, or use a plurality of the unsaturated scattered ray pixel value relationships satisfying the imaging condition in a case of imaging the subject H in combination. In addition, it is preferable to calculate the unsaturated scattered ray pixel value by, for example, “unsaturated scattered ray pixel value=unsaturated pixel value obtained in a case in which the pixel value is not saturated×scattered ray fraction in the direct radiation region”.

As described above, the region detection unit 150 distinguishes the subject region 160 and the direct radiation region 161 by using the region threshold value, but the direct radiation region 161 detected by the region threshold value may include the region transmitted through the soft region close to the skin in the subject H. In a case in which the scattered ray is also emitted from such a soft region, in a case in which the scattered ray image is subjected to the subtraction processing from the first and second radiation images G1 and G2, the scattered ray component is excessively removed. Therefore, a boundary position between the subject region 160 and the direct radiation region 161 can be adjusted by a specific width, and the scattered ray image can be obtained based on the region detection image and the scattered ray spread information after adjusting the boundary position.

The boundary position adjustment unit 154 adjusts the boundary position between the subject region 160 and the direct radiation region 161 by the specific width. It is preferable to perform adjustment of the boundary position in a case in which the pixel value of the direct radiation region exceeds the pixel value threshold value or in a case in which the dose of the radiation exceeds a dose threshold value. In addition, it is preferable to determine the specific width based on at least one of the part of the subject H, the imaging method for the subject H, or the imaging conditions. In addition, it is preferable to determine the specific width based on the pixel value in the vicinity of the boundary position in the radiation image.

Specifically, in a case in which the specific width is determined from the shape of the subject H among the parts of the subject H, that is, the body thickness distribution of the subject H, as shown in the body thickness distribution of the specific line 162 portion in FIG. 52, a portion Px in the direct radiation region 161 detected by the region threshold value, in which the body thickness exceeds 0, is used as the specific width. By adjusting a boundary position BP by the specific width corresponding to the portion Px, the subject region 160 is expanded by the specific width, while the direct radiation region 161 is narrowed by the specific width. In addition, in a case in which the specific width is determined based on the pixel value in the vicinity of the boundary position, it is preferable to determine the specific width based on the pixel-to-pixel distance between the pixel at the boundary position and the pixel within a specific range from the boundary position, which is a vicinity pixel of which the pixel value has the specific range.

In the tenth embodiment, the subtraction unit 63 derives the bone part image Gb and the soft part image Gs by using the processed first and second radiation images G1 and G2.

Note that the scattered ray removal processing according to the tenth embodiment may be performed in the first to fifth and eighth to ninth embodiments.

In addition, in each of the embodiments described above, the results of estimation of the bone part image Gb and the soft part image are derived from the simple radiation image G0, but the present disclosure is not limited to this. For example, the results of estimation of the muscle image and the fat image may be derived. Hereinafter, this case will be described as an eleventh embodiment.

FIG. 53 is a diagram showing a functional configuration of an information derivation device according to the eleventh embodiment. Note that in FIG. 53, the same reference numerals are assigned to the same configurations as those in FIG. 7, and the detailed description thereof will be omitted. In the eleventh embodiment of the present disclosure, the muscle image Gm and the fat image Gf are derived instead of deriving the bone part image Gb and the soft part image Gs as the correct answer data 42. Therefore, as shown in FIG. 53, an information derivation device 50G according to the eleventh embodiment further comprises a muscle image derivation unit 110 and a fat image derivation unit 111 with respect to the information derivation device 50 according to the first embodiment.

The muscle image derivation unit 110 derives the muscle image Gm by using the soft part image Gs. Therefore, the muscle image derivation unit 110 derives the muscle mass for each pixel in the soft region in the soft part image Gs based on the pixel value. The soft tissue includes the muscle tissue, the fat tissue, the blood, and the water. In the muscle image derivation unit 110 according to the eleventh embodiment, a tissue other than the fat tissue in the soft tissue is regarded as the muscle tissue. That is, in the muscle image derivation unit 110 according to the eleventh embodiment, a non-fat tissue including the blood and the water in the muscle tissue is handled as the muscle tissue.

The muscle image derivation unit 110 separates the muscle and the fat from the soft part image Gs by using a difference in an energy characteristic between the muscle tissue and the fat tissue. As shown in FIG. 54, the dose of the radiation after being transmitted through the subject H is lower than the dose of the radiation before being incident on the subject H, which is a human body. In addition, since the energy absorbed by the muscle tissue and the energy absorbed by the fat tissue is different and attenuation coefficients are different, the energy spectra of the radiation after being transmitted through the muscle tissue and the radiation after being transmitted through the fat tissue in the radiation after being transmitted through the subject H are different. As shown in FIG. 54, the energy spectrum of the radiation transmitted through the subject H and emitted to each of the first radiation detector 5 and the second radiation detector 6 depends on a body composition of the subject H, specifically, a ratio between the muscle tissue and the fat tissue. Since the fat tissue is more likely to transmit the radiation than the muscle tissue, the dose of the radiation after being transmitted through the human body is smaller in a case in which the ratio of the muscle tissue is larger than the ratio of the fat tissue.

Therefore, the muscle image derivation unit 110 separates the muscle and the fat from the soft part image Gs by using the difference in the energy characteristic between the muscle tissue and the fat tissue described above. That is, the muscle image derivation unit 110 generates a muscle image and a fat image from the soft part image Gs.

Note that a specific method by which the muscle image derivation unit 110 separates the muscle and the fat from the soft part image Gs is not limited, but as an example, the muscle image derivation unit 110 according to the eleventh embodiment derives the muscle image from the soft part image Gs by Expression (28) and Expression (29). Specifically, first, the muscle image derivation unit 110 derives a muscle ratio rm(x,y) at each pixel position (x,y) in the soft part image Gs by Expression (28). Note that in Expression (28), μm is a weighting coefficient depending on an attenuation coefficient of the muscle tissue, and μf is a weighting coefficient depending on an attenuation coefficient of the fat tissue. In addition, Δ(x,y) indicates a concentration difference distribution. The concentration difference distribution is a distribution of a concentration change on the image, which is seen from a concentration obtained by making the radiation reach the first radiation detector 5 and the second radiation detector 6 without transmitted through the subject H. The distribution of the concentration change on the image is calculated by subtracting the concentration of each pixel in the region of the subject H from the concentration in a blank region obtained by directly emitting the radiation in the soft part image Gs to the first radiation detector 5 and the second radiation detector 6.

rm(x,y)={μf−Δ(x,y)/T(x,y)}/(μf−μm)  (28)

Moreover, the muscle image derivation unit 110 generates a muscle image Gm from the soft part image Gs by Expression (29). Note that x and y in Expression (29) are the coordinates of each pixel of the muscle image Gm.

Gm(x,y)=rm(x,y)×Gs(x,y)  (29)

The fat image derivation unit 111 generates the fat image Gf from the soft part image Gs by using the muscle ratio rm by Expression (30).

Gf(x,y)=(1−rm(x,y))×Gs(x,y)  (30)

In the eleventh embodiment, the muscle image Gm and the fat image Gf, which are derived by the information derivation device 50Q are used as the correct answer data of the teacher data. FIG. 55 is a diagram showing the teacher data derived in the eleventh embodiment. As shown in FIG. 55, the teacher data 40A consists of the learning data 41 including the first and second radiation images G1 and G2, and correct answer data 42A including the muscle image Gm and the fat image Gf.

By learning the neural network by using the teacher data 40A shown in FIG. 55, it is possible to construct the learned neural network 23A that outputs the muscle image Gm and the fat image Gf as the results of estimation in a case in which the simple radiation image G0 is input.

Note that in the first to tenth embodiments, the results of estimation of the bone part image Gb and the soft part image Gs are derived, but the present disclosure is not limited to this. The result of estimation of any one of the bone part image Gb or the soft part image Gs may be derived. In this case, the learned neural network 23A need only be constructed by performing learning with the teacher data in which the correct answer data is any one of the bone part image Gb or the soft part image Gs.

In addition, in the eleventh embodiment, the results of estimation of the muscle image Gm and the fat image Gf are derived, but the present disclosure is not limited to this. The result of estimation of any one of the muscle image Gm or the fat image Gf may be derived. In this case, the learned neural network 23A need only be constructed by performing learning with the teacher data in which the correct answer data is any one of the muscle image Gm or the fat image Gf.

Note that in each of the embodiments described above, the estimation device 10 learns the neural network to construct the learned neural network 23A, but the present disclosure is not limited to this. The learned neural network 23A constructed in a device other than the estimation device 10 may be used for the estimation unit 23 of the estimation device 10 in the present embodiment.

In addition, in each of the embodiments described above, the first and second radiation images G1 and G2 are acquired by the one-shot method in a case in which the energy subtraction processing is performed, but the present disclosure is not limited to this. The first and second radiation images G1 and G2 may be acquired by a so-called two-shot method in which imaging is performed twice by using only one radiation detector. In a case of the two-shot method, a position of the subject H included in the first radiation image G1 and the second radiation image G2 may shift due to a body movement of the subject H. Therefore, in the first radiation image G1 and the second radiation image G2, it is preferable to perform the processing according to the present embodiment after registration of the subject is performed.

As registration processing, for example, a method disclosed in JP2011-255060A can be used. In the method disclosed in JP2011-255060A, for each of the first and second radiation images G1 and G2, a plurality of first band images and a plurality of second band images representing structures having different frequency bands are generated, a misregistration amount of the positions corresponding to each other in the first band image and the second band image of the corresponding frequency band is acquired, and the registration of the first radiation image G1 and the second radiation image G2 is performed based on the misregistration amount.

In addition, in each of the embodiments described above, the derivation of the bone part image Gb and the soft part image Gs as the correct answer data of the teacher data and the estimation processing of the bone part image Gb and the soft part image Gs from the simple radiation image G0 are performed by using the radiation image acquired by the system that images the subject H by using the first and second radiation detectors 5 and 6, it is needless to say that the technology of the present disclosure can be applied to even in a case in which the first and second radiation images G1 and G2 are acquired by using an accumulative phosphor sheet instead of the radiation detector. In this case, the first and second radiation images G1 and G2 need only be acquired by stacking two accumulative phosphor sheets, emitting the radiation transmitted through the subject H, accumulating and recording radiation image information of the subject H in each of the accumulative phosphor sheets, and photoelectrically reading the radiation image information from each of the accumulative phosphor sheets. Note that the two-shot method may also be used in a case in which the first and second radiation images G1 and G2 are acquired by using the accumulative phosphor sheet.

In addition, the radiation in the embodiments described above is not particularly limited, and α-rays or γ-rays can be used in addition to X-rays.

In addition, in the embodiments described above, various processors shown below can be used as the hardware structures of processing units that execute various pieces of processing, such as the image acquisition unit 21, the information acquisition unit 22, the estimation unit 23, the learning unit 24, and the display controller 25 of the estimation device 10, and the image acquisition unit 61, the scattered ray removal unit 62, and the subtraction unit 63 of the information derivation device 50. As described above, the various processors include, in addition to the CPU that is a general-purpose processor which executes software (program) and functions as various processing units, a programmable logic device (PLD) that is a processor whose circuit configuration can be changed after manufacture, such as a field programmable gate array (FPGA), and a dedicated electric circuit that is a processor having a circuit configuration which is designed for exclusive use in order to execute specific processing, such as an application specific integrated circuit (ASIC).

One processing unit may be configured by one of these various processors, or may be a combination of two or more processors of the same type or different types (for example, a combination of a plurality of FPGAs or a combination of the CPU and the FPGA). In addition, a plurality of the processing units may be configured by one processor.

As an example of configuring the plurality of processing units by one processor, first, as represented by a computer, such as a client and a server, there is an aspect in which one processor is configured by a combination of one or more CPUs and software and this processor functions as a plurality of processing units. Second, as represented by a system on chip (SoC) or the like, there is an aspect of using a processor that realizes the function of the entire system including the plurality of processing units by one integrated circuit (IC) chip. In this way, as the hardware structure, the various processing units are configured by using one or more of the various processors described above.

Moreover, as the hardware structures of these various processors, more specifically, it is possible to use an electrical circuit (circuitry) in which circuit elements, such as semiconductor elements, are combined. 

What is claimed is:
 1. An estimation device comprising: at least one processor, wherein the processor functions as a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple two-dimensional image acquired by simply imaging the subject, and the learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images.
 2. The estimation device according to claim 1, wherein the emphasis image for learning is derived by energy subtraction processing of performing weighting subtraction on the two radiation images.
 3. The estimation device according to claim 2, wherein the emphasis image is at least one of a bone part image in which a bone part of the subject is emphasized or a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning is derived by recognizing the bone part and the soft part of the subject by using at least one radiation image of the two radiation images, deriving attenuation coefficients of the bone part and the soft part by using results of recognition of the bone part and the soft part and the two radiation images, and performing the energy subtraction processing by using the attenuation coefficients.
 4. The estimation device according to claim 2, wherein the emphasis image is a bone part image in which a bone part of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning is derived by deriving new weighting coefficients used for the weighting subtraction based on a pixel value of the bone part included in the bone part image and a pixel value of the soft part included in the soft part image, deriving a new bone part image and a new soft part image by performing the weighting subtraction on the two radiation images by using the new weighting coefficients, and repeating derivation of a further new weighting coefficient based on the new bone part image, derivation of a further new weighting coefficient based on the new soft part image, and derivation of a further new bone part image and a further new soft part image based on the further new weighting coefficients.
 5. The estimation device according to claim 2, wherein the emphasis image is a bone part image in which a bone part of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, and the emphasis image for learning is derived by deriving, for each of different energy distributions, a difference between a value of an attenuation coefficient of the soft part×a thickness of the soft part+an attenuation coefficient of the bone part×a thickness of the bone part, and each pixel value of the radiation image while changing, from initial values, the attenuation coefficient of the soft part for each of different energy distributions, the thickness of the soft part, the attenuation coefficient of the bone part for each of different energy distributions, and the thickness of the bone part, deriving the attenuation coefficient of the soft part and the attenuation coefficient of the bone part for each of different energy distributions, at which the difference is minimized or the difference is smaller than a predetermined threshold value, and performing the energy subtraction processing by using a weighting coefficient derived based on the attenuation coefficient of the soft part and the attenuation coefficient of the bone part.
 6. The estimation device according to claim 2, wherein the emphasis image for learning is derived by deriving a composition ratio of a plurality of compositions included in a soft part of the subject, deriving, for each pixel of the two radiation images, an attenuation coefficient for each of different energy distributions depending on the composition ratio, and performing the energy subtraction processing by using a weighting coefficient derived based on the derived attenuation coefficient.
 7. The estimation device according to claim 6, wherein the composition ratio is obtained by deriving a body thickness of the subject as a first body thickness and a second body thickness for each pixel of each of the two radiation images and deriving the composition ratio for each pixel of the radiation image based on the first body thickness and the second body thickness.
 8. The estimation device according to claim 7, wherein the composition ratio is obtained by deriving the first body thickness and the second body thickness based on an attenuation coefficient of each of the plurality of compositions for each of different energy distributions, deriving the first body thickness and the second body thickness while changing a thickness of the composition and the attenuation coefficient of each composition, and deriving the composition ratio based on the thickness of the composition in which a difference between the first body thickness and the second body thickness is equal to or smaller than a predetermined threshold value.
 9. The estimation device according to claim 2, wherein the emphasis image for learning is derived by performing scattered ray removal processing of removing a scattered ray component of the radiation emitted to the subject, which is scattered by the subject, from the two radiation images, and performing the energy subtraction processing on the two radiation images from which the scattered ray component is removed.
 10. The estimation device according to claim 9, wherein the scattered ray removal processing is performed by acquiring a radiation characteristic of an object interposed between the subject and a radiation detector that detects the radiation image depending on the body thickness distribution of the subject, deriving a primary ray distribution and a scattered ray distribution of the radiation included in each of the two radiation images by using the imaging condition, the body thickness distribution, and the radiation characteristic of the object, deriving an error between a sum of the primary ray distribution and the scattered ray distribution of each of the two radiation images and a pixel value at each position of the two radiation images, updating the body thickness distribution such that the error is smaller than a predetermined threshold value, repeating derivation of the radiation characteristic based on the updated body thickness distribution and derivation of the primary ray distribution and the scattered ray distribution included in each of the two radiation images, and subtracting the scattered ray distribution in a case in which the error is smaller than the predetermined threshold value from each of the two radiation images.
 11. The estimation device according to claim 9, wherein the scattered ray removal processing is performed by deriving a first primary ray distribution and a scattered ray distribution of the radiation transmitted through the subject by using the two radiation images, deriving a second primary ray distribution and a scattered ray distribution of the radiation transmitted through an object interposed between the subject and a radiation detector that detects the radiation image by using the first primary ray distribution, the scattered ray distribution, and a radiation characteristic of the object, and deriving the radiation images after transmission through the subject and the object by using the second primary ray distribution and the scattered ray distribution.
 12. The estimation device according to claim 9, wherein the scattered ray removal processing is performed by deriving a region detection image by detecting a subject region in which the radiation is transmitted through the subject and reaches a radiation detector and a direct radiation region in which the radiation directly reaches the radiation detector without being transmitted through the subject in the two radiation images, deriving a scattered ray image relating to the scattered ray component based on the region detection image and scattered ray spread information relating to spread of a scattered ray, and subtracting the scattered ray image from the two radiation images.
 13. The estimation device according to claim 2, wherein the emphasis image for learning is obtained by deriving a processing content of first granulation suppression processing on a first radiation image having higher S/N among the two radiation images, deriving a processing content of second granulation suppression processing on a second radiation image having lower S/N based on the processing content of the first granulation suppression processing, performing granulation suppression processing on the first radiation image based on the processing content of the first granulation suppression processing, performing granulation suppression processing on the second radiation image based on the processing content of the second granulation suppression processing, and deriving the emphasis image for learning by using the two radiation images subjected to the granulation suppression processing.
 14. The estimation device according to claim 13, wherein the processing content of the first granulation suppression processing is derived based on a physical quantity map of the subject derived based on at least one of the first radiation image or the second radiation image.
 15. An estimation method comprising: using a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple radiation image acquired by simply imaging the subject to derive the result of estimation of the at least one emphasis image in which the specific composition of the subject is emphasized from the simple radiation image, wherein the learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images.
 16. A non-transitory computer-readable storage medium that stores an estimation program causing a computer to execute a procedure of: using a learned neural network that derives a result of estimation of at least one emphasis image in which a specific composition of a subject including a plurality of compositions is emphasized from a simple radiation image acquired by simply imaging the subject to derive the result of estimation of the at least one emphasis image in which the specific composition of the subject is emphasized from the simple radiation image, wherein the learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject with radiation having different energy distributions and an emphasis image for learning in which the specific composition of the subject is emphasized, which is derived from the two radiation images. 